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FOREWARD FOR USER'S GUIDE-V 


The Flexible Spacecraft Dynamics Program (FSD program) is a second generation computer 
program having evolved from the Radio Astronomy Explorer (RAE) Dynamics program 
which had its inception in 1965. The main FSD program was written and developed by the 
AVCO Corporation, Systems Division, Wilmington, Massachusetts and was dciivoicd ‘ 

in 1970. Since that time many checkouts and numerous improvements have been made. The 
program has run on the IBM 360/91 and 360/95 machines at GSFC, requiring up to 600 
bytes of storage with overlay capability. More recently, the program has been running on the 
IBM 3081 machine which has no restrictions on storage requirements. Currently, the use of the 
program has been extended to the VAX 11/780 machine. The operation of the program is user 
oriented. That is. the dcMgn of the program and inpur output is such that, aside from initial job 
control language cards, a dynamics! or an analytical control engineer can set up and run pro- 
blems without programming assistance: no programming skills are required. 

It is beyond the scope of this Guide to go into the many applications of this program. Exami- 
nation of the input control words and various options does give an indication of its versatility. 

The program can be used in dynamics and control analysis as well as in orbit support of 
deployment and control of spacecraft. This program has been used to simulate the dynamics 
of antenna deployment and in-orbit attitude performance of the RAE-A. B. IMP-I. J, ISEE-A. 

C. Langley HAWKEYE. Air Force SCATHA. Italy SAN MARCO-D. Japan EXOS-B. ISPM. 
Dynamic Explorer-A. B. ISTP and CRRES. 

With current emphasis on active control of pointed instruments on flexible spacecraft, 
recent additions have been in the controls area. Some additions to the current FSD program 
are as follows: 

1. A two axis gimble platform and digital control system (proportional, integrai and 
derivative controller) to track the earth’s magnetic field from a spinning or non-spin- 
ning spacecraft. The magnetic field can be sensed on the spacecraft hub or at the end 
of a flexible boom. The control system coding is contained in one subroutine so that 
other control laws can be implemented without requiring changes in the overall 
program. 

2. A two axis platform magnetic tracking system using the pm-uni n i.pmimmi d.smpci 
degree of freedom for the azimuth gimble. This permits simulation of flexible ele- 
ments on the azimuth platform as well as on the spacecraft while tracking the mag- 
netic field. 

3. Proportional, integral and derivative controller has been added to each axis » n v : pisj* 
and yaw of the earth pointing mode) to generate control torques from the respective 
momentum wheels. An arbitrarily oriented momentum wheel with control is also 
included. User formulated control laws are permitted. 

4. Jet damping added to the thruster option to simulate launch vehicle dynamics. 

5. Thermal expansion and contraction of flexible elements caused by changing sun 
angle, spacecraft or orbit shadowing. 
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In addition to being used in Goddard related flight programs, the FSD program 
is being modified and maintain ed by the Systems Division of the AVCO Corpora- 
tion, Wilmington, Massachusetts. 

J. V. Fedor 
E. A. Law lor 1 
J. F. Downey* 

A. H. Forbes : 
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^Old Dominion Systems. Inc. 
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A USER’S GUIDE TO THE FLEXIBLE SPACECRAFT 
DYNAMICS AND CONTROL PROGRAM-V 


INTRODUCTION 

This computer program was initially developed to simulate the dynamics of the 
IMP and RAE class of satellites. Generality was retained in its development 
so that it is applicable to the simulation of the dynamics and control of a large 
class of flexible and rigid spacecraft. 

Hie program is applicable to inertially oriented spinning, earth oriented or 
gravity gradient stabilized spacecraft. Internal and external environmental ef- 
fects developed at orbital altitude are simulated. The effects include gravity 
gradient forces, solar pressure, magnetic torques and thermal bending due to 
solar heating. Body torquing devices in the computer program include momen- 
tum wheels, a viscous ring nutation damper, magnetic torquer coils and attitude 
control thrusters. For gravity gradient satellites, an option is available for 
simulating either a magnetic hysteresis or viscous libration damper. 

The computer program has the capability of simulating up to ten flexible tubular 
elements arbitrarily oriented with respect to the body fixed coordinate frame. 

A finite series of shape function are used to describe the bending and twisting of the 
flexible elements. Higher order displacement terms are retained in order to achieve 
reasonable accuracy for large displacements. 

The equations of motion are derived from variational principles, i.e. , the prin- 
ciple of virtual work. The generalized coordinates include the three rotational 
and three translational degrees of freedom of the body fixed axes and the ampli- 
tudes of the shape functions for each flexible element. An additional generalized 
coordinate is necessary to describe the morions of the libration damper. 

Generalized forces were derived and programmed for gravity gradient forces , 
solar pressure, bending stresses and structural damping. The induced tempera- 
ture gradients and solar pressure generalized forces are derived from the in- 
stantaneous angle of incidence between the sun line vector and the deformed 
flexible elements. The effects of aerodynamic drag on the flexible motions of 
a spacecraft is also computed for low altitude orbits. 

A special purpose computer program (Integral Evaluation Program) was also 
developed to compute input data for the dynamics computer program. This com- 
puter program evaluates definite integrals that evolve in the mathematical proc- 
ess of spatially integrating the internal and external forces acting on the flexible 
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elements of a satellite. The integrals are normalized products of the shape 
functions and their derivatives evaluated over the flexible elements lengths. 

For a given shape function, selected to represent the deformed shape of a flexi- 
ble element, the integrals have to be evaluated only once. The integrals are 
read into the dynamics program either on cards or compiled into block data. 

The dynamics simulations can then be made without further recourse to the in- 
tegral evaluation program. The shape functions are specified by the coefficients 
of polynomials. For a flexible element with no tip mass, a set of typical shape 
function used would represent cantilever beam bending modes. Other more ap- 
propriate shape functions would be specified for simulating flexible elements 
with tip masses. The dynamics program can use up to three shape functions or 
modes in simulation of the deflectional motion of the flexible elements. At times 
it is necessary to have different types of flexible elements with different stiff- 
ness characteristics on the same spacecraft. The dynamics program has, 
therefore, the capability of utilizing two different sets of values as determined 
by the integral evaluation program for two different families of shape function. 
An example of the use of this capability would be a spinning spacecraft requiring 
interlocked closed cross-section elements on the spin axis and utilizing wire 
elements on the transverse axes. 

The input to the Flexible Spacecraft Dynamics program consists of four main 
parts, i.e.: 

1. Input which is necessary to construct the orbit of the spacecraft. 

2. Control words to invoke or delete various options such as gravity gra- 
dient effects, control torques , starting integration interval, etc. 

3. Input to describe mathematically the spacecraft, appendages and control system. 

4. Desired state variable output which is controlled by 1NOPT 1 or 2 and by 
KPLOTS array addresses. 

The input to these portions is on punched cards , punched in columns 2 through 
72. The Fortran input symbol can be punched in any of these columns followed 
by its input value. Since the program utilizes a nameless read, more than one 
input symbol and associated value can be punched on a card provided the symbol 
value pairs are separated by at least one blank. Within each main portion, the 
order of input is immaterial. A heading card can be inserted within the input 
deck provided an H is punched in column one. Comment cards can be inserted 
in the input dat eck with an "$" or "*" punched in column one. 

The input to the main parts will be described in detail. 
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PART 1 


ORBIT GENERATION 


There are two options for describing the orbit. They are: 

1. Utilizing the internal orbit generator which generates a two-body 
(Keplerian) orbit at each time step in the program. Constant orbit 
drifting rates of longitude of ascending node and argument of perigee 
can also be input to the program providing better orbit simulation if 
long term study is required. For the majority of simulations, this 
simplified orbit is quite adequate. However, if a definitive orbit is 
desired, or the program is used during operational support of a mis- 
sion, the second orbit option should be used. 

2. An orbit tape can be read by the program via control words and Fortran 
I/O unit. The orbit tape is generated by GTDS (Goddard Trajectory 
Determination System) of Code 582 at GSFC for the satellite and time 
span in question. 

The orbit option is set by two input control words, IORB and 1TAPE. 

For internal orbit generation, set 

IOR3 = 0, ITAPE = 0 (preset values) 

For reading orbit tape option, set 

IORB = 1, ITAPE = NM 

where NM is the Fortran I/O unit number corresponding to the fol- 
lowing JCL to be inserted for this option 

//GO. FTNMF001 DD UNIT=2400,DSN=FSD. EPHEM, 
//VOL-=SER=XXXXXX, LABEL=1, BLP), DCB= 

(RECFM=VS, BLKSIZE=2808, BUFNO=l ) 

In either case the following orbital input parameters are input: the parameters 
for which the earth is the central body are pre-set internally, so it is not really 
necessary to read in these values unless it is desired to change them. The user 
can also specify a central body other than the earth by reading in the appropriate 
parameters. 


3 




The orbital input parameters are as follows: 


EARTH PROPERTIES 


Fortran 

Symbol 

Math 

Symbol 

Description 

Units 

AEARTH 

«E 

Semi -major axis of Earth. Internally 
set to 6378. 165 

km 

ZMU 

% 

Gravitational constant of Earth. In- 
ternally set to 3.986U32D + 05 

km 3 /sec 2 

ZJ20 


Oblate Earth coefficient. Internally 
set to 1.082D-06 

N.D. 

FLAT 

a 

a-b 

a = Semi-major axis of Earth 
b = Semi-minor axis of Earth 
Flat is internally set to 298. 3 

N.D. 

wwo 


GHA of Aries at instant vernal equinox 
occurs in the year of simulation. Used 
for the (2, 2) term in the expansion of 
Earth's gravitational model (preset = 55.0) 

WE 


Earth's mean angular velocity about 
the sun. Internally set to 1. 1407D- 
05 Deg/ sec . 

deg/ sec 

TVER 


Time in seconds from Jan. 1st at 
00:00 to time of Vernal Equinox in the 
year ot simulation. Internally set to 
6873720. 0 seconds. 

sec 

ECLPTC 


Obliquity ->f the ecliptic. Internally set 
to 23. 444 deg. 


J2 

J3 


Oblate Earth coefficients. Internally 
set as follows: 


J4 


J2 = 1082. 3D-06, J3 = 2.3D-06, 
J4 = 1.8D-06, J22 = 5.35D-06 


J22 


Non-dimensional values. 
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Classical Orbit Parameter Option 


Fortran 

Symbol 

Math 

Symbol 

Description 


Units 

IKPLR 

N. A. 

Control word for Keplerian orbit 
parameter input 

N. D. 



IKPLR = 0, input | Xj | and |x, f 
IKPLR = 1, (preset) input Keplerian 
orbit parameters 


Fortran 

Symbol 

Math 

Symbol 

Description 

Preset 

Value 

Units 

AS 

a 

Semi-major axis 

8000.0 

km 

E 

e 

Eccentricity 

0.0 

N.D. 

El 

i 

Inclination angle 

0.0 

deg 

F 

f 

True anomaly 

0.0 

deg 

BW 

ft 

Longitude of ascending node 

0.0 

deg 

W 

CO 

Argument of perigee 

0.0 

deg 

BWDOT 

ft 

Rate of change of longitude 
of ascending node 

0.0 

deg/day 

WDOT 

CO 

Rate of change of argument 
of perigee 

0.0 

deg/day 

The following cards are input only if IKPLR = 0, otherwise they are 

omitted. 

XSAT (1-3) 


Components of initial position 


km 


vector in the Aries (Equatorial) 
Inertial Frame. 

(preset = 8000.0, 0.0, 0.0) 


XSATDT (1-3) i {x,} Components of initial velocity km/sec 

vector in the Aries (Equatorial) 

Inertial Frame. 

(preset - 0.0, 7.0587, 0.0) 
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PART 2 


CONTROL WORDS 

TIME CONTROL WORDS 
Fortran 

Symbol Description Units 

IDATE Year, Month, Day YYMMDD (This number should fall Integers 

within the orbit described in Part 1). (i.e. it should 
be the same year as the value of TVER specified, 
preset = 760101) 

TIME Time in seconds from start of day. This is also the Seconds 

problem start time (preset = 0.0) 

TSTOP Problem stop time. Specified as incremental time Seconds 

from TIME (preset = 3600.0) 

FREQ Output data print frequency (preset = 60.0) Seconds 

GENERAL CONTROL WORDS 
Fortran 

Symbol Description Units 

IDATA Control word for printing of input data by sub- 
routine ECHO (preset = 1) 

If IDATA = 0 prints out original set of inpul parameters, implies 
MLAST = 1. 

If IDATA = 1 prints out set of input parameters which were the end 
conditions of the previous case in the stack. Implies 
MLAST = 0. 

IDATA, unlike MLAST, is used in the case within the stack to which 
it applies. 

MLAST Control word associated with stacking of Integer, 

cases. If a card MLAST = 0 is inserted N.D. 

within a case , the terminal conditions of this 
case are used as initial conditions for the 
next case which follows. 
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Fortran 

Sjymbol Description Units 

If a card MLAST = 1 is inserted in each case, N.D. 

the first set of data for spacecraft descrip- 
tion is used in each case which follows the 
first except for those cards which override 
the originals iu the first case. (Preset = 0). 

NOTL 

Stacking Procedure 

The capability is provided to stack computer runs in order to simplify 
and speed the parametric analysis of spacecraft dynamics problem y. 

Two types of problems can be processed: (1) Parallel stacking capa- 
bility — series of simulations are being made with similar initial con- 
ditions and it is desired to change one or more parameters in each case, 
or (2) Series stacking capability — a continuous (in time) simulation is 
run where at specified times, changes in parameters are made but the 
dependent variables are carried over from the last time of one se- 
quence to the beginning of the next sequence. 

To run a stacked case, it is only necessary to insert a "1" card, i.e. , 
a "1" punched in column one, after each set of input data. The simu- 
lation will continue until the TSTOP of a sequence is reached. The 
new data will be read in -fter the "1" card and the simulation will re- 
start with the modified input data. An indefinite number of stacked 
cases can be processed providing there is sufficient computer time 
available. 

. deployment run is also considered to be a stacked case. The first 
sequences should have either MDPLY = 1 or DDPLY = 1 or both. De- 
ployment is either initiated or terminated in a particular sequence by 
setting ZLl(k) equal to the desired deployment rate or equal to zero, 
respectively. 

Note that it is important to specify an appropriate value of DELTAT 
for each sequence in order to start the Runge Kutta-Adams Moulton 
integration scheme. If no new DELTAT is specified in the new stack, 
the DELTAT value of the previous stack will be used lor the new 
stack. 

If it is desired to rerun the deployment sequence with a change in a 
parameter, set MLAST = 1 in the last sequence. TMs will restore 
the input set of initial conditions. 


MLAST 

(con’t) 
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OPTION CONTROL WORDS 


Fortran 

Symbol 

INOPT 

IHAMLT 

IGRAV 

MDPLY 

DDPLY 


ISDPLY 


Description 

Control word for type of satellite being studied. 

(Preset value = 1) 

INOPT = 1 Spin Stabilized Satellite 

INOPT = 2 Gravity-Gradient or earth-pointing Sat- llite 

Control word to call Relative Hamiltonian subroutine. 
(Preset value = 0) 

IHAMLT = 0 By-passes subroutine 

IHAMLT = 1 Employs subroutine 

Control word to incorporate or exclude gravity 
effects. (Preset value = 1) 

IGRAV = 0 Gravity excluded 

IGRAV = 1 Gravity included 

Control word for boom deployment subroutine. 

(Preset vaiue = 0) 

MDPLY = 0 Not a deployment case 

MDPLY = 1 Deployment of booms required 

Control word for libration damper boom 
deployment. (Preset value = 0) 

DDPLY = 0 Not a deployment case 

DDPLY = 1 Deployment of libration damper boom 
required 

Control word to begin boom deployment from a "see 
sun" pulse. (Preset value = 0) 

ISDPLY = 0 By-passes subroutine 

ISDPLY = 1 Employs subroutine 


Units 

Integer, 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 
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Fortran 

Symbol 

IDAMP 

mscs 

ISPIN3 

IATTDE 

IWHEEL 

IMGMTS 

ITORK 


Description 


Control word to exclude or allow libration damper 
motion. (Preset value = 0) 

IDAMP = 0 No damper motion 

IDAMP = 1 Damper motion 

Control word to employ viscous ring nutation 
damper into program. (Preset value = 0) 

IVISCS = 0 By-passes subroutine 

IVISCS = I Employs subroutine 

Control word to call spin axis moment subroutine. 
(Preset value = 0) 

ISPIN3 = 0 Bypass subroutine 

ISPIN3 = 1 Employs subroutine 

Control word to call the attitude control subroutines. 
(Preset value ■= 0) 

IATTDE -- 0 Bypasses subroutine 

IATTDE = 1 Employs subroutine 

Control word to call momentum wheel subroutine. 
(Preset value = 0) 

IWHEEL - 0 Bypasses subroutine 

IWHEEL = 1 Employs subroutine 

Control word to exclude or allow magnetic moments. 
(Preset value = 0) 

IMGMTS = 0 magnetic moments excluded 
IMGMTS = 1 magnetic moments allowed 

Control word to apply torque about body axes. 
(Preset value = 0) 

ITORK = 0 No torque applied 

ITORK = 1 Torque applied 


Units 

Integer, 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 
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Fortran 

Symbol 

IBENDM 
ITENSE 
DP LANS 


BATSH 


IWRTTF 

IAFM(l) 

IAFM(2) 


Description 


Control word to calculate and print out flexible ele- 
ment root bending moments. (Preset value = 0) 

IBENDM = 0 No calculation 

IBENDM = 1 Calculate and print 

Control word to calculate and print out element 
root tensions. (Preset value = 0) 

ITENSE = 0 No calculation 

ITENSE = 1 Calculate and print 

Control word to invoke a thermal lag effect on 
appendages due to planet and satellite shadowing. 
(Preset value = 0) 

IP LANS = 0 No lag invoked 

IP LANS = 1 Lag invoked 

Control word to invoke thermal lag effect of 
appendages due to satellite shadowing 
(Preset value = 0) 

ISATSH = 0 No shadowing 

ISATS1I = 1 Shadowing and lag 

Control word to print out thermal forces before and 
after computation of thermal lag. (Preset value = 0) 

IWRTTF = 0 No print out 

IWRTTF « l Print 

Control word to print out unit sun vector in body frame. 
(Preset value = 0) 

IAFM(l) = 0 No sun vector printout 

LAFM(l) = 1 FTints out sun vector 

Control word for punch cards restart (Preset 
value = 0) 


Units 

Integer, 

N.D. 


Integer, 

N.D. 


Integer , 
N.D. 


Integer, 

N.D. 


Integer, 

N.D. 


Integer, 
N. D. 


Integer, 

N.D. 
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Fortran 

Symbol 


IAFM(3) 


IAFM(4) 


IAFM(5) 


IACOMP 


fflUBAC 


ITIPAC 


Description 


IAFM(2) = 0 No cards punched 

IAFM(2) = 1 Cards punched 

Control word for printing of output data. 

(Preset value = 0) 

IAFM(3) = I No pr—tout, printer plots only 

IAFM(3) = 0 Printout and plots 

Control word to activate the data records and auto- 
correlation function plots for fast fourier transform 
analysis. (Preset value = 0) 

IAFM(4) = 0 No data record and autocorrelation 
function plots 

IAFM(4) = 1 Data record and autocorrelation 
function plots 

Control word to activate Adams -Moultion integrator 
message table. (Preset value = 0) 

IAFM(5) = 0 No integrator message 

IAFM(5) = 1 Integrator message printed 

Control word to activate the acceleration computation 
for both hub and element tip accelerations. 

(Preset value = 0) 

IACOMP = 0 Bypasses computation 

IACOMP = 1 Computes accelerations 

Control word to activate acceleration computation 
and printout for the hub. (Preset value = 0) 

IHUBAC = 0 No computation 

IHUBAC - 1 Computes hub accelerations 

Control word to activate acceleration computation 
and printout for the element tip. (Preset value = 0) 

ITIPAC = 0 No computation 

ITIPAC = 1 Computes tip accelerations 


Units 


integer, 

N.D. 


Integer, 

N.D. 


Integer, 
N. D. 


-iteger, 

N.D. 


Integer, 
N. D. 


Integer, 

N.D. 
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Fortran 

Symbol 

IPVLSE 

ISPLSE 

IPLPRP 

ISPNP 

IHCALC 


IHREF 


IKPLR 


Description 

Control word to activate thrust subroutine. 

(Preset value * 0) 

IPULSE = 0 Bypasses subroutine 

IPULSE > 1 Employs subroutine 

Control word to activate sun crossing time to start 
thrusting. (Preset value = 0) 

ISPLSE = 0 Sun crossing not used 

ISPLSE - 1 Sun crossing used 

Control word for number of thrust pulses per spin period. 
IPLPRP = 1 One pulse (Preset) 

IPLPRP = 2 Two pulses 

Control word to print out orbit update message (only 
if ISPLSE = 1, IPULSE > 1). i.e., if ISPNP = 5, (preset) 
the orbit update message will be printed at every 5th 
pulse. 

Control word to calculate and print out the angular 
momentum vector of the system and system 
moments of inertia (preset value = 0) 

IHCALC = 0 No calculation 

IHCALC = 1 Calculates and prints 

Control word to compute reference direction for 
EPSERR from the initial system angular 
momentum vector (Preset value = 0) 

IHREF = 0 No calculation 

IHREF = 1 Determines reference direction from 
initial angular momentum vector 

Control word for Keplerian orbit parameter input 
(Preset value = 1) 

IKPLR = 0 Input position an velocity 
IKPLR = 1 Input Keplerian orbit parameters 


Units 

Integer, 
N. D. 


Integer, 
N. D. 


Integer. 

N.D. 


Integer, 
N. D. 


Integer, 
N. D. 


Integer, 

N.D. 


Integer, 
N. D. 
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Fortran 

Symbol 

ICSD 

IPRY 

ISCP 

1PUNCH 

KNTRL(IO) 

I2BDY 

IRAST 

1ARSTO) 


Description 

Control word to activate the Fast Fourier Transform analysis 
and plot. See p. 18 for KPLOTS control also (Preset value = 0) 

ICSD = 0 By-passes subroutine 

ICSD = 1 Employs subroutine 

Control word to print and plot pitch, roll and yaw rates instead 
of OMBC values when INOPT = 2 

IPRY = 0 No print out and plot 

IPRY = 1 Print out and plots (preset = 0) 

Control word to activate the spacecraft configuration plot op- 
tion. (Not available in current program) 

ISCP = 0 By-passes subroutine 

ISCP = 1 Employs plot subroutine 

Control word for punched card output under normal case end. 
Set MLAST = 1 and IPUNCH = I (preset = 0) 

Vector of control integers ( 10) for DE-B control system 
simulation 

Control word to add secondary body to simulation 

I2BDY = 0 No secondary body 

I2BDY = 1 Secondary body present 

Control word to invoke prescribed rastering motions for 
secondary body. 

IRAST = 0 No rastering 

IRAST = 1 Rastering Prescribed 

Control word to specify type of rastering cycle to be invoked 
on each axis. 

IARST (I) = 0 No motion 
IARST (I) = 1 Motion of Type 1 
IARST (I) = 2 Motion of Type 2 
I = 1 to 3 

I = 1 motion about the 3 axis 


Units 

Integer. 

N.D. 


Integer. 

N.D. 


Integer. 

N.D. 

Integer. 

N.D. 

Integer. 

N.D. 


Integer. 

N.D. 


Integer, 

N.D. 
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Fortran 

Symbol 


Description 


Units 


IACFLT(20) 

IDPHLD 

IPLOT 

IGMBL 

IPCTCS (20) 
IJTOMP 

IGMBLD 

IPLDCS (20) 
ITHRM (20) 
lAMPRM(l) 


I = 2 motion about the 1 axis 
I = 3 motion about the 2 axis 

Control integers for filtered sensor signal to actuate momentum 
wheel cycling or pulsed thrusting. 

Control word for simulation of constant angular velocity of 
libration damper. 

IDPHLD = 0 No simulation 

IDPHLD = I Constant angular velocity for damper 
Control word for plotting of output data. 

IPLOT = 0 No plots 

IPLOT - 1 Printer plots or FFT analysis. (Preset = 0) 

Control word for two axis gimble simulation. 

IGMBL 0 No gimble simulated (preset) 

IGMBL 1 Gimble simulated 

20 Control words for gimble platform con rol system 
simulation (preset = 0) 

Control word to invoke jet damping. 

IJTDMP 0 No jet damping (preset) 

IJTDMP 1 Jet damping 

Control word for two axis gimble (damper simulation) 

IGMBLD 0 No gimble simulated (preset) 

IGMBLD 1 Gimble simulated 

20 Control words for gimble platform (damper) simulation 
(preset = 0) 

20 Control words for element thermal expansion - contraction 
simulation 

Control word for arbitrary oriented momentum wheel 
simulation (preset = 0) 
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Integer, 

N.D. 

Integer. 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 

Integer, 

N.D. 


Integer, 

ND. 

Integer, 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 


Integer, 

N.D. 



Fortran 

Symbol 


Description 


Units 


1AMWH(10) 

10 Control words for axis momentum wheel simulation 
(preset = 0) 

Integer, 

N.D. 

lOMKDM (1) 

Control word to invoke use of reference frequency for damping 
coefficient (preset = 0) 

Integer, 

N.D. 

INTEGRATION CONTROL WORDS 


Fortran 

Symbol 

Description 

Units 

DELTAT 

Starting integration interval for Runge Kutta. 
(preset value = 0.01) 

Seconds 

DELMIT 

Minimum value for time integration interval for 
predictor-corrector. If integration step size be- 
comes less than DELMIT, simulation stops. 
(Preset = 1. 0 D-7) 

Seconds 

FACTOR 

Percentage by which At is varied in the Adams- 
Moulton integrator. If the difference between extrapo- 
lated and interpolated values of the integrands is 
greater than the upper bound (on the difference) At is 
decreased; if the difference is greater than the lower 
bound, A t is increased, (predictor-corrector) 
FACTOR is presently set to 0.3). i.e., At n+J = 

At n (1 * FACTOR) 

N.D. 
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Description Units 

The upper bound on the absolute difference between 
the extrapolated and interpolated values used in the 
predictor-corrector of the integration routine. 

The lower bound on the absolute difference between 
the extrapolated and interpolated values. 

NOTE 

UP and DN are bounds on the difference between the extrapolated and inter- 
polated values of the components of the state vector. They are set internally 
and are not required to be input. Unless the user fully understands their use, 
it is recommended that they be left undisturbed. The same applies to the 
integration bounds control words which follow. 

INTEGRATION BOUNDS 


Corresponds 
to specific 
variables 

Corresponds 
to specific 
variables 


Fortran 

Symbol 

UP(i) 

DN(i) 


Fortran 

Symbol Description Units 


CONSA 


COMEG 


DUC 


This is the upper bound on the absolute difference N. D. 

between the extrapolated and interpolated values 
used in the predictor-corrector for the integra- 
tion of the elements in the first two rows of di- 
rection cosine matrix. The lower bound is auto- 
matically set two orders of magnitude below the 
upper bound. (Preset value * 1. 0 D-5) 

This is the upper bound allowed for the error in rad/sec 

the predictor-corrector oi the integration of the 
body components of inertial angular velocity. 

The lower bound is automatically set two orders 
of magnitude below the upper bound. 

(Preset value = 1. 0 D-7) 

This is the upper bound on the error for the rad 

predictor-corrector integration for the libration 
damper angle of gravity-gradient satellites. For 
the simulation of spin stabilized satellites with 
nutation damping, this represents the upper 
bound on the error for predictor-corrector inte- 
gration for the viscous torque of the damper. 

(Preset value = i. 0 D-2) 
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Fortran 

Symbol 

Description 

Units 

DUC1 

This is a lower bound for DUC. 
(Preset value = 1.0 D-4) 

rad 

DUCD 

This is the upper bound on the error for the 
predictor-corrector integration for the angular 
velocity of the libration damper for gravity - 
gradient satellites. For the simulation of spin- 
stabilized satellites with viscous damping, this is 
the upper bound on the error for the predictor- 
corrector integration for the angular momentum 
of the viscous liquid. (Preset value = 1.0 D-3) 

rad/ sec 

DUCD1 

This is the lower bound for DUCD. 
(Preset value = 1.0 D-5) 

rad/sec 

DOOP 

This is the upper bound on the error for predictor- 
corrector integration for the 2 axis components 
of damper element deflections. 

(Preset value = 1.0 D-l) 

feet 

DOOP1 

This is the lower bound for DOOP. 
(Preset value = 1.0 D-3) 

feet 

DOOPV 

This is the upper bound on the error for predictor- 
corrector integration for the 2 axis components of 
damper element velocity. (Preset value = 5.0 D-5) 

ft/ sec 

DOOPV1 

This is the lower bound for DOOPV. 
(Preset value = 5.0 D-7) 

ft/sec 

DIP 

This is the upper bound on the error for predictor- 
corrector integration for the 3 axis components of 
damper element deflections. 

(Preset value = 1.0 D-l) 

feet 

DIP1 

This is the lower bound for DIP. 
(Preset value = 1.0 D-3) 

feet 
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Fortran 

Symbol 

DIPV 

DIPV1 

AOOP 

AOOP1 

AOOPV 

AOOPV1 

AIP 

AIP1 

AIPV 

AIPV1 

TWIUP 


Description 


This is the upper bound on the error for predictor- 
corrector integration for the 3 axis components of 
damper element velocities. 

(Preset value - 5.0 D-7) 

This is the lower bound for DIPV. 

This is the upper bound on the error for the 
predictor-corrector integration for the 2 axis 
element frame component of antenna tip 
displacement. (Preset value - 1.0 D-l) 

This is the lower bound for AOOP. 

(Preset value - 1.0 D-3) 

This is the upper b.oynd on the error for the 
predictor-corrector integration for the 2 axis 
element frame component of tip velocities. 

(Preset value - 5.0 D-5) 

This is the lower bound for AOOPV. 

(Preset value = 5.0 D-7' 

This is the upper bound on the error for the 
predictor-corrector integration for the 3 axis 
element frame component of tip displacement. 
(Preset value - 1.0 D-l) 

This is the lower bound for AIP. 

(Preset value - 1.0 D-3) 

This is the upper bound on the error for the 
predictor-corrector integration for the 3 axis 
element frame component of tip velocity. 

(Preset value - 5.0 D-5) 

This is the lower bound for AIPV. 

(Preset value - 5.0 D-7) 

Upper integration bound for twist angle. 

(Preset value - 1. 0 D-4) 


Units 

ft/sec 


ft/sec 

feet 


feet 

ft/sec 


ft/sec 

feet 


feet 


ft/sec 

deg 
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Fortran 

Symbol 

TWIDN 

TWDUP 

TWDDN 

CSUP(20) 

CSDN(20) 

SBUP(2) 

SBDN(2) 

ACPARM( 1 9) 
ACPARM(20) 


Description 


Lower integration bound for twist angle. 

(Preset value = 1.0 D-6) 

Upper integration bound for twist velocity. 

(Preset value = 1.0 D-5) 

Lower integration bound for twist velocity 
(Preset value = 1.0 D-7) 

Upper bound on difference between predicted and corrected 
control system state vector. Location in CSUP corresponds 
to the location of the variable in the state vector initial 
condition array SVCS. 

(Preset value = 1.0 D-02) 

Lower bound on difference between predicted and corrected 
confol system state vector. Location in CSDN corresponds 
to location of the variable in the state vector initial 
condition array SVCS. 

(Preset value = 1.0 D-04) 

Integration upper bounds for secondary body angles and 
angular rates. 

(Preset = 1.0 D-3) 

Integration lower bounds for secondary body angles and 
angular rates. 

(Preset = 1.0 D-5) 

Upper bound for filter integrater. 

Lower bound for filter integrater. 


Units 

deg 

deg/sec 

deg/sec 


rad 

rad/sec 

rad 

rad/sec 
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Fortran 

Symbol 

Description 

Units 

PCSPRM(30) 

Integration upper bound for magnetometer first order lag 
transfer function (preset = 0) 

gauss 

PCSPRM(31) 

Integration lower bound for magnetometer first order lag 
transfer function (preset = 0) 

gauss 

GMUP(1) 

Integration upper bound for gimble angles (preset = 1.0 D-3) 

rad 

GMUP(2) 

Integration upper bound for gimble angular rates - i i) 3i 

rad/sec 

GMDN(l) 

Integration lower bound for gimble angles (preset = 1.0 D-5) 

rad 

GMDN(2) 

Integration lower bound for gimble angular rates (preset = 1.0 D-5) 

rad/sec 

DMUP(l) 

Integration upper bound for damper gimble angles (preset = 1.0 D-3) 

rad 

DMUP(2) 

Integration upper bound for damper gimble anguiar rates 
(preset = 1.0 D-3) 

rad/sec 

DMDN(l) 

Integration lower bound for damper gimble angles (preset = 1.0 D-5) 

rad 

DMDN(2) 

Integration lower bound for damper gimble angular rates 
(preset = 1.0 D-5) 

rad/sec 

DCSPRM(30) 

Integration upper bound for damper related magnetometer first 
order lag transfer function (preset = 0) 

gauss 

DCSPRM(31) 

Integration lower bound for damper related magnetometer first 
order lag transfer function (preset = 0) 

gauss 

THRMPR(2) 

Integration upper bound for thermal expansion - contraction 
simulation (preset = 0) 

°R 

THRMPR (3) 

Integration lower bound tor thermal expansion-contraction simulation (preset 

= 0) °F 

AMPARM 

(101) 

Integration upper bound for filtered error angle of the arbitrarily 
oriented momentum wheel (preset = 0) 

deg 

AMPARM 

(102) 

Integration lower bound for filtered error angle of the arbitrarily 
oriented momentum wheel (preset = uj 

deg 

AMPARM 

(103) 

Integration upper bound for arbitrarily oriented momentum 
wheel speed (preset — 0) 

deg/ser 

AMPARM 

(104) 

Integration lower bound for arbitrarily oriented momentum 
wheel speed (preset = 0) 

deg/ sec 

AMWHPR 

(101) 

Integration upper bound for filtered angle error of the axis 
momentum wheel simulation (preset = 0) 

deg 
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AMWHPR 

(102) 

Integration lower bound for filtered angle error of the axis 
momentum wheel simulation (preset = 0) 

4 -* 

AMWHPR 

(103) 

Integration upper bound for axis momentum wheel speed 
(preset = 0) 

dcv.'sec 

AMWHPR 

(104) 

Integration lower bound for axis momentum wheel speed 
(preset = 0) 

deg/sec 

PLOT CONTROL WORDS 


Fortran 

Symbol 

Description 

Units 

IPLOT 

Control word for plotting of output data. 

Integer, 

N.D. 


IPLOT = 0 No plots. 



IPLOT = 1 Printer plots or FFT analysis. 
(Preset = 0). 

Integer, 

N.D. 
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Fortran 

Symbol Description 


Preset 

Value Units 


IPLMOD Control word for individual modal analysis 0 Integer 

plot N.D. 

IPLMOD = 0, No individual modal 
analysis plot 

IPLMOD = 1, Store individual modal 
variables for plotting 


SCMOD Control word to activate the flexible element 1 Integer 

requiring independent higher mode plots. E.g., N.D. 

if NELMTS = 7, NDAMPR = 0, Set IPLMOD = 1 
and IKMOD = 5 in order to activate higher 
modal displacement plots for 5th elements 


NOTE 

Value of MODES for IKMODth element must be 
greater or equal to one and less or equal to 3. 


The variables to be plotted by printer plot are controlled by KPLOTS array 
values. The dimensions of KPLOTS array is 253. The control of plotting is as 
follows 


KPLOTS(I) = 0 
KPLOTS(I) = 1 
KPLOTS(I) = 2 

KPLOTS(I) = 3 


do not plot Ith variable (preset value) 

plot Ith variable 

for plot and FFT analysis 

for FFT analysis only 


All KPLOTS array addresses are integer variables with value of either 0 or 1. 
For the first ten KPLOTS array addresses, (i.e. , 1 = 1 to 10), the variables 
plotted are functions of INOPT. For I 3* 11, the plotting variables are inde- 
pendent of the INOPT value. The KPLOTS array addresses are given in the 
following table. 
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For I < 10 


KPLOTS Array 

Fortran Variable Plotted vs. Time 


Address 

INOPT 1 

INOPT 2 

Units 

1 

PSI1 

ALFAE 

deg 

2 

THET1 

BETAE 

deg 

3 

PKI1 

GAMAE 

deg 

4 

OMEG1 

OMBC1 

deg/sec 

5 

OMEG2 

OMBC2 

deg/ sec 

6 

OMEG3 

OMBC3 

deg/sec 

7 

PHILD 

PHILD 

deg 

8 

NUT ANG (deg) 

PRAT 

deg/ sec 

9 

EPSERR (deg) 

RRAT 

deg/ sec 

10 

— 

YRAT 

deg/ sec 


The definition of the Fortran variables are given in Part 4 of this document. 
For I > 11 


K PLOTS 
Array 

Address Fortran 


I 

Description 

Symboi 

Units 

Note 

11 

1st flexible element 1-2 
plane tip displacement 

U21 

feet 


12 

1st flexible element 1-2 
plane tip velocity 

U2DOT1 

ft/ sec 


13 

1st flexible element 1-3 
plane tip displacement 

U31 

feet 
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KPLOTS 

Array 

Address 

I Description 


14 

1st flexible element 1-3 
plane tip velocity 

U3DOT1 

ft/ sec 

15 

1st flexible element twist* 
displacement 

no 

symbol 

deg 

16 

1st flexible element twist* 
velocity 

no 

symbol 

deg/ sec 

17 

2nd flexible element 1-2 
plane tip displacement 

U22 

feet 

18 

2nd flexible element 1-2 
plane tip velocity 

U2DOT2 

ft/sec 

19 

2nd flexible element 1-3 
plane tip displacement 

U32 

feet 

20 

• 

2nd flexible element 1-3 
plane tip velocity 

U3DOT2 

ft/sec 

• 

* 

• 

ft 

etc. 



65 

10th flexible element 1-2 
plane tip displacement 

U210 

feet 

66 

10th flexible element 1-2 
plane tip velocity 

U2DOT10 

ft/sec 

67 

10th flexible element 1-3 
plane tip displacement 

U310 

feet 

68 

10th flexible element 1-3 
plane tip velocity 

U3DOT10 

ft/sec 

69 

10th flexible element* 
twist displacement 

no 

symbol 

deg 


‘lndniclu.il Ivusl modes are printed out and the sunt ot the modes is plotted. See p. I()h lor output print sunbeds 


Fortran 

Symbol Units Note 
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KPLOTS 

Array 

Address 

1 

Description 

Fortran 

Symbol 

Units 

Note 

70 

10th flexible element* 
twist velocity 

no symbol 

deg/sec 


71 

1st damper flexible ele- 
ment 1-2 plane tip 
displacement 

UD21 

feet 


72 

1st damper flexible ele- 
ment 1-2 plane tip velocity 

UD2DT1 

ft/sec 


73 

1st damper flexible ele- 
ment 1-3 plane tip 
displacement 

UD31 

feet 


74 

1st damper flexible ele- 
ment 1-3 plane tip velocity 

UD3DT1 

ft/sec 


75 

1st damper element* 
twist displacement 

no symbol 

deg 


76 

1st damper element* 
twist velocity 

no symbol 

deg/sec 


77 

2nd damper flexible ele- 
ment 1-2 plane tip 
displacement 

UD22 

feet 



etc. 




128 

10th damper flexible ele- 
ment 1-3 plan j Jp velocity 

UD3DTJ0 

ft/sec 


129 

10th damper element* 
twist displacement 

no symbol 

deg 


130 

10th damper element* 
twist velocity 

no symbol 

deg/sec 


131 

Mode 1 1-2 plane tip dis- 
placement for element 
IKMOD 

AK1 

feet 

1 < MODES < 3 
IPLMOD = 1 


‘Individual twist modes arc primed out and the sunt of the modes is plotted. See p. 106 for output print symbols. 
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K PLOTS 
Array 
Address 
I 

Description 

Fortran 

Symbol 

Units 

Note 

132 

Mode 1 1-3 plane tip 
displacement for 
element IKMOD 

BK1 

feet 

1 < MODES < 3 
IPLMOD = 1 

133 

Mode 2 1-2 plane tip 
displacement for 
element IKMOD 

AK2 

feet 

1 < MODES < 3 
IPLMOD = 1 

134 

Mode 2 1-3 plane tip 
displacement for 
element IKMOD 

BK2 

feet 

1 < MODES < 3 
IPLMOD = 1 

135 

Mode 3 1-2 plane tip 
displacement for 
element IKMOD 

AK3 

feet 

1 < MODES < 3 
IPLMOD = 1 

136 

Mode 3 1-3 plane tip 
displacement for 
element IKMOD 

BK3 

feet 

1 < MODES < 3 
IPLMOD = 1 

137 

Mode 1 twist displace- 
ment for element IKMOD 


deg 


137 

Mode 1 twist displace- 
ment for element IKMOD 

CWK1 

deg 


138 

Mode 2 twist displace- 
ment for element IKMOD 

CWK2 

deg 


139 

Mode 3 twist displace- 
ment for element IKMOD 

CWK3 

deg 


140 

Mode 1 1-2 plane tip dis- 
placement for damper 
element IKMOD 


feet 


141 

Mode 1 1-3 plane tip dis- 
placement for damper 
element IKMOD 

DINK2 

feet 


142 

Mode 2 1-2 plane tip dis- 
placement for damper 
element IKMOD 

DINK 2 

feet 
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K PLOTS 
Array 
Address 

I Description 

Fortran 

Symbol 

Units 

Note 

143 

Mode 2 1-3 plane tip dis- 
placement for damper 
element IKMOD 

DINK2 

feet 


144 

Mode 3 1-2 plane tip dis- 
placement for damper 
element IKMOD 

DOUTK3 

feet 


145 

Mode 3 1-3 plane tip dis- 
placement for damper 
element IKMOD 

DINK3 

feet 


146 

Mode 1 twist displace- 
ment for damper element 
ZKMOD 

CWK1 

deg 


147 

Mode 2 twist displace- 
ment for damper element 
ZKMOD 

CWK2 

deg 


148 

Mode 3 twist displace- 
ment for damper element 
ZKMOD 

CWK3 

deg 


149 

154, 

Accelerometer reading 
value from 1 to 6 
► respectively 

ACCRED1 

ACCRED2 

ACCRED6 

ft/sec 2 

ft/sec 2 

• 

• 

• 

ft/ sec 2 

IACOMP = 1 
IHUBAC = 1 
The Max. No. 
of plot is 
NUMHUB 

155 

Instantaneous body 
moment of inertia about 
1-axis 

BDCX 

slug-ft 2 

IHCALC = 1 

156 

Instantaneous body 
moment of inertia about 
2 -axis 

BIYY 

slug-ft 2 

IHCALC = 1 

157 

Instantaneous body 
moment of inertia about 
3-axis 

BIZZ 

slug-ft 2 

IHCALC = 1 

158 

Right ascension angle of 
angular momentum vector 
in Aries inertial frame 

None 

27 

deg 

IHCALC = 1 
PLOTTED ONLY 



K PLOTS 
Array 
Address 
I 


159 


160 

161 



171 

172 

173 

174 
» 

189 

190 

191 

192 


Fortran 

Description Symbol Note 


Declination angle of 
angular momentum vector 
in Aries inertial frame 

None 

deg 

IHCALC * 1 
PLOTTED ONLY 

Magnitude of angular 
momentum vector 

HMAG 

ft-lb-sec 

IHCALC = 1 

Root tension history 

TENSN1 

lb 

ITENSE = 1 

for element 1 to 10 

TENSN2 

lb 

plotting quantities 

respectively 

TEN3N10 

• 

< 

(NELMTS+ NDAMPR) 

1-2 plane root bending 
moment for element 1 

BNMTA1 

lb- ft 

IBENDM = 1 

1-3 plane root bending 
moment for element 1 

BNMTB1 

lb-ft 

EBENDM = 1 

1-2 plane root bending 
moment for element 2 

BNMTA2 

lb-ft 

IBENDM = 1 

1-3 plane root bending 
moment for element 2 

etc. 

BNMTB2 

lb-ft 

IBENDM = 1 

1-2 plane root bending 
moment for element 10 

BNMTA10 

lb-ft 

IBENDM = 1 

1-3 plane root bending 
moment for elemen, 10 

BNMTB10 

lb-ft 

IBENDM = 1 

Right ascension angle of 
the unit sun line vector 
in body reference frame 

None 

deg 

IAFM(l) = 1 
PLOTTED ONLY 

Declination angle of the 
unit sun line vector in 
body reference frame 

None 

deg 

IAFM(l) - 1 
PLOTTED ONLY 
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KPLOTS 

Array 

Address 

I 

193 ^ 

> 

195 > 

196 

197 

206 

207 

208 

209 

210 
211 
212 


Description 

Fortran 

Symbol 

Units 

Note 


Magnetic field vector 
components in 1, 2, 3, 
body axis respectively 

SMAGB1 

SMAGB2 

SMAGB3 

GAUSS 

GAUSS 

GAUSS 

IMGMTS » 

1 

Relative Hamiltonian of 
the system 

HAMILT 

ft-lb 

IHAMLT = 

1 

Magnitude of root bending 
moment for 1st flexible 
element 

None 

ft- lb 

PLOTTED ONLY 

etc. 





Magnitude of root bending 
moment for 10th flexible 
element 

None 

ft-lb 

PLOTTED ONLY 

Momentum wheel speed 
1-axis 

WHL SPD1 

deg/ sec 



Momentum wheel speed 
2-axis 

WHL SPD2 deg/ sec 



Momentum wheel speed 
3-axis 

WHL SPD3 deg/sec 



External moment about 
1 body axis 

MOMENT 1 

ft-lbs 



External moment about 
2 body axis 

MOMENT 2 

ft-lbs 



External moment about 
3 body axis 

MOMENT 3 

ft-lbs 
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KPLOTS 




Array 




Address 

Fortran 



I Description 

Symbol 

Units 

Note 


213 Component of angular HBODY 1 ft-lb-sec IHCALC 1 
momentum on 1 body axis 


214 

Component of angular 
momentum on 2 body axis 

HBODY 2 

ft-lb-sec 

IHCALC 1 

215 

Component of angular 
momentum of 3 body axis 

HBODY 3 

ft-lb-sec 

IHCALC 1 

216 

Pitch sensor output 


volts 

IWHEEL 1 
and 

217-219 

Pitch sensor dynamics 



KNTRL(l) 
1 or 2 

220 

Not used 




221 

Roll sensor output 


volts 


222-224 

Roll sensor dynamics 


volts/sec, 

etc 


225 

Not used 




226 

Output of pitch compen- 
sation amplifier 


volts 


227 

Not used 




228 

Tachometer output 


volts 


229 

Not used 




230 

Momentum wheel speed 




231-233 

Not used 




234 

Nutation damper phase shift 
dynamics 




235 

Nutation damper phase shift 





output 
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KPLOTS 

Array 

Address 

I 

Description 

Fortran 

Symbol 

Units 

236 

Rotation of secondary 
body about 3 axis 

GAMSB 

deg 

237 

Rotation of secondary 
body about the carried 1 
axis 

ALPSB 

deg 

238 

Rotation of secondary 
body about the carried 2 
axis 

BETSB 

deg 

2*9 

Relative angular rate of 
secondary body about 3 
axis 

GAMMAD 

deg/sec 

240 

Relative angular rates of 
secondary body about the 
carried 1 axis 

ALPHAD 

deg/sec 

241 

Relative angular rates of 
secondary body about the 
carried 2 axis 

BETAD 

deg/sec 

242 

Component of secondary 
body relative angular velo- 
city on 1 axis of secondary 
body 

OM1SB 

deg/sec 

243 

Component of secondary 
body relative angular velo- 
city on 2 axis of secondary 
body 

OM2SB 

deg/sec 

244 

Component of secondary 
body relative angular velo- 
city on 3 axis of secondary 
body 

OM3SB 

deg/sec 

245 

State variable sensor output 


varies 

246 

State variable filter output 


varies 

247 

State variable filter first inte 
grater output 

31 

varies 


Note 



KPLOTS 

Array 

Address 

I 

Description 

Fortran 

Symbol 

Units 

Note 

248 

Product of inertia 

IXY 

slug ft 2 


249 

Product of inertia 

IXZ 

slug ft 2 


250 

Product of inertia 

IXY 

slug ft 2 


251 

1st damper element 
1-2 plane slope 


rad 


252 

1 st damper element 
1-3 plane slope 


rad 


253 

• 

2nd damper element 
1-2 plane slope 


rad 


• 

• 

etc. 


rad 


269 

10th damper element 
1-2 plane slope 


rad 


270 

10th damper element 
1-3 plane slope 


rad 


271 

1st flexible element 
1-2 plane slope 

U2P 1 

rad 


272 

1st flexible element 
1-3 plane slope 

U3P2 

rad 


273 

2nd flexible element 
1-2 plane slope 

U2P 2 

rad 



etc. 




289 

10th flexible element 
1-2 plane slope 

U2P 10 

rad 


290 

10th flexible element 
1-3 plane slope 

U3P 10 

rad 


291 

292 

1st mode 1-2 plane slope 
1st mode 1-3 plane slope 

32 

rad 

rad 

Mode output is for 
the element specified 
by input IKMOD 



KPLOTS 

Array 


Address 

I 

Fortran 

Description Symbol 

Units 

Note 

293 

2nd mode 1-2 plane slope 

rad 


294 

2nd mode 1-3 plane slope 

rad 


295 

3rd mode 1-2 plane slope 

rad 


296 

3rd mode 1-3 plane slope 

rad 


297 

Azimuth angle for the gimble 
azimuth platform 

deg 


298 

Azimuth angular rate for the 
gimble azimuth platform 

deg/sec 


299 

Elevation angle for the gimble 
elevation platform 

deg 


300 

Elevation angular rate for the 
gimble elevation platform 

deg/sec 


301 

Output of first order lag transfer 
function for magnetometer one 
axis 

Gauss 


302 

Output of first order lag transfer 
function for magnetometer two 
axis 

Gauss 


303 

Output of first order lag transfer 
function for magnetometer 
three axis 

Gauss 


304 

Azimuth error output 

deg 

IPLTCS 1 

305 

Elevation error output 

deg 

IPLTCS 1 


306 Azimuth PID digital controller 
output 

Units depend on PID gain 
constant units 

307 Elevation PID digital controller 
output 



33 



KPLOTS 

Array 

Address 

I 

Fortran 

Description Symbol 

Units 

Not 

308 

Azimuth platform drive motor 
torque 

ft-lbs 

IPLTCS 1 

309 

Elevation plailnnn drive mc*or 
torque 

ft-lbs 


310 

Azimuth angle for the gimble 
azimuth platform 

deg 


311 

Azimuth angular rate for the 
azimuth platform 

deg/sec 

IGMBLD 1 

312 

Elevation angle for the gimble 
elevatiOii platform 

deg 


313 

Elevation angular rate for the 
gimble elevation platform 

deg/sec 


314 

Output of first order lag transfer 
function for magnetometer one 
axis 

gauss 


315 

Output of first order lag transfer 
function for magnetometer two 
axis. 

gauss 


316 

Output of first order lag transfer 
function for magnetometer 
three axis 

gauss 


317 

Azimuth error output 

rad 


318 

Elevation error output 

rad 


319 

Azimuth PID digital controller 
output 


Units depend on 
PID gain constants 

320 

1. luxation PID digital controllor output 



321 

322 

Azimuth platform drive motor 
torque 

Elevation platform drive motor 
torque 

ft-lbs 

ft-lbs 

IAMWH(1) = 1 
and 

IAMWH(2) = 1 
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KPLOTS 




Array 




Address 

Fortran 



I Description 

Symbol 

Units 

Note 


323 Roll axis filter output of error rad 

signal for momentum wheel 
control 


324 

Pitch axis filter output of error 
signal for momentum wheel 
control 

rad 

IAMWH(3) = 1 

325 

Yaw axis filter output of error 
signal for momentum wheel 
control 

rad 

IAMWH(4) - 1 

326 

Roll axis momentum wheel speed 

deg/sec 


327 

Pitch axis momentum wheel speed 

deg/sec 


328 

Yaw axis momentum wheel speed 

deg/sec 


329 

Arbitrarily oriented momentum wheel 
filter output of control system 
angular error 


lAyoPMd)-- 1 

330 

Momentum wheel speed 



331 

1st element temperature 

°F 

ITHRM(l) = 1 
and 

ITHRMfl 1) = 1 


etc 


etc 

340 

10th element temperature 

°F 

1THRM(20) = 1 


35 


Fort i an 

Symbol Description 

LCPU CPU time in seconds allowed for print out, integer 

(Preset = 20). 

DIAGNOSTIC CONTROL WORDS 

Fortran Preset 

Symbol Description Value 

IOUT Control word to print computations as performed in the 1 

Simulator. 

IOUT = 1 No print of computations. 

IOUT = 2 Print of computations at print frequency. 

IOUT = 3 Print of computations for each call of 
DEREQ1. 

KLUGE Control word to stop computing on this input set. If 0 

KLUGE = 1 and IOUT = 2, machine stops at time zero, 
before normal output is printed, and prints 
calculations. 

ISWTCH Control word to print energy related computations. If 8 

ISWITCH - 0 and IHAMLT = 1, energy computations 
are printed out. If ISWITCH ^ 0, no print of energy 
computations are made. 

NOPT A control word to give the number of times subroutine 0 

DEREQ1 was entered. 

NOPT = 0 No printout of number of times. 

NOPT = 1 Prints out number of times . 


Present 

Value 
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PART 3 

SPACECRAFT DESCRIPTION & SIMULATION OPTIONS 
CORE PROPERTIES 


Fortran 

Symbol Description Units 


BDYMI (i, j) Iy Moments of inertia of satellite core about slug ft 2 

Y body axes, i = 1, 2, 3 j = 1, 2, 3 

SCO S 0 Projected area of central core of satellite ft 2 

Used for aerodynamics and solar pressure. 

(preset = 14.6) 

ZMS Mass of entire satellite (preset = 25.0) slugs 

HUBCDA (3) CpA values for spacecraft hub along three ft 2 

body axes, (preset = 1. 0, 1. 0, 1. 0) 

HUBCP (3) Body frame position vector of hub center of ft 2 

pressure, (preset = 0. 0, 0. 0, 0. 0) 


ELEMENT GEOMETRY AND PHYSICAL PROPERTIES 


Fortran 

Symbol 

Math 

Sjymbol 

Description 

Units 

NELMTS 


The number of elements rigidly attached 
to satellite core, (preset = 6, do not set to 0). 

Integer, 

N.D. 

ZL0 (k) 

®0k 

Length of the k* element at start of 
problem time, (preset = 1.0, do not set to 0) 

feet 

ZL1 (k) 


Velocity of deployment of k 111 element, 
(preset = 0.0) 

ft/sec 

ZLA (k) 

*k 

Acceleration of deployment of k 1 * 1 
element, (preset = 0.0) 

ft/sec 2 

ZBZ (l,k)2 
ZBZ (2, k) 2 
ZBZ (3,k)2 

^lk 

z m2k 

z m3k 

Coordinates of the origin of the k ele- 
ment frames as defined ir. the refer- 
ence Y body frame or iteration damper 
Z frame, (preset - 0.0) 

feet 
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Symbol 


Math 

Symbol 


Description 


Units 


ALFAEK(k) <* k 


BETAEK(k) 0 k 


GAMAEK(k) 




MODES (k) 


A (k , j) 1 A~ 


ADOT (k, j) 1 Ajj 


B(k,j)l Bjj 


Euler angle relating orientation of k 1 * 1 deg 

element frame as defined in the refer- 
ence Y body frame or libration damper 
Z frame. (Around 1 axis, 2nd angle 
in 2 - 1 - 3) rotation, (preset = 0.0) 

Euler angle relating orientation of k’ n deg 

element frame as defined in the ref- 
erence Y body frame or libration 
damper Z frame. (Around 2 axis, 1st 
angle in 2 - 1 - 3 rotation.) (preset = 0.0) 

Euler angle relating orientation of the deg 

k* element frame as defined in the 
reference Y body frame or libration 
damper Z frame. (Around 3 axis, 

3rd angle in 2 - 1 - 3 rotation.) (preset - 0.0) 

Designation of bending modes for Integers 

element 

0 = Rigid body 

1 = Bending mode 1 

2 = Bending mode 2 

3 = Bending mode 3 (includes 3 modes) 

Component of the k * core element tip feet 

deflection in the j * bending mode as 
measured along the 2 axis of the ele- 
ment frame, (preset = 0.0) 

Component of the k** 1 core element tip ft/sec 
velocity in the j 1 * 1 bending mode as 
measured in the element frame along 
the 2 axis of the frame, (preset = 0.0) 

Component of the k 1 * 1 core element tip feet 

deflection in the j * bending mode as 
measured in the element frame along 
the 3 axis of the frame, (preset = 0.0) 
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Math 

Symbol Description Units 

B- Component ol the k ,h core element tip ft/sec 

velocity in the bending mode as 
measured in the element frame along 
the 3 axis of the frame, (preset = 0.0) 

NOTE 

In the above arrays, the "1" following the parentheses indi- 
cates that the first integer in the parentheses varies as num- 
bers are read off across the input card. For example, the 
card: 


Fortran 

Symbol 

BDOT (k, j) 1 


A(k, 1)1 2.0 4.0 5.0 

would be interpreted by the program as: 

A (1, 1) =2.0 
A (2, 1) = 4.0 
A (3, 1) = 5.0 etc. 

The arrays A, ADOT, B, BDOT are both input and output of the 
program. Thus if the user wished to start a problem on case 
with initial tip deflections and on velocities of the rigidly at- 
tached elements, these quantities would be input. Otherwise, 
they are internally set to zero at the start of the problem and 
are output only. 


Fortran 

Symbol Description 

LK (k) This is a control word for the selection of one of two preset 

(Block 2) data sets for the k 1 * 1 element. Each data set contains 
normalized mass integrals and normalized effective areas inte- 
grals for flexible elements. The first data set is generated 
using cantilever beam modes. The second data set utilizes 
spinning string modes with a tip mass. 

LK (k) = 1 uses data set 1 

LK (k) = 2 uses data set 2 
(preset to 1) 
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Symbol Description 


A data set can be generated for a specific antenna configuration 
by using the Integral Evaluation Program. 

LLK (k) This is a control word for the selection of one of two preset data 

sets for the k* element. Each data set contains normalized in- 
ternal force integrals and normalized thermal force integrals for 
flexible elements . 

LLK (k) = 1 uses data set 1 


Fortran 

LLK (k) = 2 uses data set 2 
(preset to 1) 

A data set can be generated for a specific flexible configuration 
by using the Integral Evaluation Program. 

Math 

Symbol 

Symbol 

Description 

Units 

EMODLS (k) 

E 

Young modulus for k th element 
(preset = 2.0 D7) 

lb- in. 2 

RTUBE <k) 

r 

Mean radius of k 111 element 
(preset = 2.935 D-l) 

inches 

HTUBE (k) 

b 

Wall thickness of k ,h element 
(preset = 2.0 D-3) 

inches 

THERMC (k) 

«T 

Thermal coefficient of expansion 
for k^ element, (preset = 8.85 D-6) 

in./in./°I 

TIPMS (k) 

m t 

Tip mass attached to k lh element 

slugs 

SAO(k) 

S co 

The projected area of a one foot 
length of element corrected for 
flow around a cylinder. Used in 
computing aerodynamic and solar 
pressure. (preset = 0.0) 

ft 2 

STMK (k) 

Src 

The projected area of a tip mass 
corrected for flow around a spheri- 
cal body. Used in computing aero- 
dynamic and solar pressure. 

ft 2 
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Symbol 

Math 

Symbol 

Description 

Units 

RHOK(k) 

P k 

Mass per unit length of the k 1 * 1 
element, (preset = 4.36 D-4) 

slugs/ft 

POO 

Po 

Solar pressure on a flatplate at 
normal incidence, (preset = 0.0) 

lb/ft 2 

DTOO 

AT 0 

Temperature differential across 
antenna at normal incidence of 
sunlight, (preset - 0.0) 

°F 

CDAMP(n, k)2 

C CT 

Damping ratio in n th bending mode 
for k* element, (preset = 0.0) 

N.D. 

SKOA(k, n)2 

^kno 

The 2 axis component in the k^ 
element frame of the offset zero 
stress position corresponding to 
the n* mode shape, (preset = 0.0) 

feet 

SKOB(k, n)2 

®kno 

The 3 axis component in the k th 
element frame of the offset zero 
stress position corresponding to 
the n* mode shape, (preset - 0.0) 

feet 

AERO 


Aerodynamic drag coefficients for 
the elements, (preset = 2.0) 

N.D. 

TDIS(k) 


Factor to account for variations in 
temperature distributions for k th 
element, (preset = 2.0) 

N.D. 


Computation of Flexible Element Root Bending Moments 


Fortran 

Math 



Symbol 

Symbol 

Description 

Units 

IBENDM 


Control word to calculate and print out flexible 
element root bending moments. 

(Preset value = 0) 

IBENDM = 0 No calculation 

IBENDM = 1 Calculate and print 

Integer, 

N.D. 
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Symbol Symbol Description Units 

ZKBM(6) x' n (0) Root curvatures of normalized shape function N.D. 

for first three modes. Two sets are per- 
missible corresponding to data sets controlled 
by LLK(K) (internal forces). 

LLK(K) = 1 ZKBM(l-3) (preset = 3.5, -22.0, 61.7) 
LLK(K) = 2 ZKBM(4-6) (preset = 3. 5, -22. 0, 61. 7) 

Computation of Element Root Tensions 


Fortran 

Symbol Description Units 

ITENSE Control word to calculate and print out element Integer , 

root tensions. (Preset value = 0) N.D. 

ITENSE = 0 No calculation 
ITENSE = 1 Calculate and print 

Thermal Lag-Element Bending 


Fortran Math 

Symbol Symbol Description Units 


IP LANS Control word to invoke a thermal lag effect Integer, 

on appendages due to planet and satellite N. D. 
shadowing. (Preset value = 0) 


IP LANS = 0 No lag invoked 
IP LANS = 1 Lag invoked 


TAUPL 


T p Characteristic delay time for computing sec 

planet shadowing switching times, (preset = 1.0) 


OCCRIT 


O cr Threshold for invoking thermal lag for N.D. 

planet shadowing 
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Symbol 

Math 

Symbol 

Description 

Units 

ISATSH 


Control word for invoking satellite 
shadowing and thermal lag for satellite 
shadowing. (Preset value = 0) 

Integer, 
N. D. 



ISATSH = 0 No shadowing 
ISATSH - 1 Shadowing and lag 


RADSH 

R s 

Radius of shading disk for satellite 
shadowing, (preset - 10.0) 

feet 

TAUK(IO) 

T k 

Characteristic delay time for each 
element, (preset = 1.0) 

sec 

OCULTK(IO) 

Ok 

Threshold for invoking thermal lag switch- 
ing for each element (satellite shadowing) 

N.D. 

IWRTTF 


Control word to write out thermal forces 
before and after computation of thermal 
lag. (Preset value = 0) 

Integer, 

N.D. 

Thermal Expansion and Contraction of Element 


Fortran 

Symbol 

Math 

Symbol 

Preset 

Description Value 

Units 

ITHRM(20) 


Control words for element thermal 0 

expansion simulation 

Integer, 

N.D. 

lTHRM(l) 


Basic control word for thermal 
expansion simulation 

ITHRMf 1 ) 0 No thermal expansion 

ITHRM(l) Thermal expansion 
simulated 


ITHRM(2) 


Internal. Not input 


ITHRM(3-10) 


Not used 


JTHRM( 11-20) 


Control words for individual elements 
according to internal numbering 
system 
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Fortran 

Symbol 


Math 

Symbol 


Description 


Preset 

Value 


Units 


THRMPR(IOO) 

THRMPR(l) 

THRMPR(2-3) 

THRMPR (4) 

THRMPR(5) 

THRMPR(6-10) 
THRMPR( 11-20) 

THRMPR(21-30) 

THRMPR(31-40) 

THRMPR(41-50) 

THRMPR(51-60) 


ITHRM(K+ 10) 0 No thermal 
expansion for 
KTH element 

ITHRM(K+10) 1 Thermal expansion 
for the KTH 
element 

Physical constants for thermal 0.0 DO 

expansion simulation 

Solar radiation constant at the 
satellite 


Upper and lower integration 
bounds 

Acceleration due to gravity at sea 
level 

Stefan Boltzmann Radiation 
constant 

Not used 

Spc;ific heat for individual element 
material according to internal 
numbering system 

Emissivity of individual element 
surface material according to 
internal numbering system 

Absorbtivity of individual element 
surface material according to 
internal numbering system 

Area of individual element 0.0 DO 

surface for thermal absorbtion 
according to internal numbering 
system 

Area of individual element surface 
for thermal emission according to 
internal numbering system 


Btu 

sec-ft 2 °R 

°R 

ft/sec 

Btu 

sec-ft 2 °R 4 

Btu 
lb °R 

N.D. 

N.D. 

ft 2 

ft 2 
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Thermal Expansion and Contraction of Element (Cont) 


fortran 
S\ mbol 


Math 
S\ mhol 


Preset 

Description \ due Units 


THRMPR(6 1-80) Used internally 

THRMPR(81-90) Element reference temperature for °R 

length calculation 

THRMPR(91-100) Initial condition for element °R 

temperature simulation 

TYPICAL INPUT FOR THERMAL EXPANSION SIMULATION 
ITHRM(l) 1 

THERMAL COEFFICIENT OF EXPANSION PER DEG R 
THERMCO-4) l.OD-3 

ELEMENTS IN SIMULATION 
ITHRM(ll) 10 10 

SOLAR RADIATION BTU / SEC. - FOOT* *2 
THRMPR(l) 0.12278 

GRAVITATIONAL ACCELERATION FOR CONVERSION ONLY 
THRMPR(4) 32.174 

STEPHAN BOLTZMANN CONSTANT BTU/SEC. - FOOT**2 - DEG R**4 
THRMPR(5) 0.476D-12 

SPECIFIC HEAT FOR ELEMENT MATERIAL BTU / POUND - DEG R 
THRMPRO 1-14) 0.21 

EMISSIVITY FOR ELEMENT SURFACE 
THRMPR(21-24) 0.2 

ABSORPTIVITY FOR ELEMENT SURFACE 
THRMPR(31-34) 0.2 

AREA ASSOCIATED WITH ABSORPTIVITY FOOT**2 
THRMPR(4 1 -44) 0.05 

AREA ASSOCIATED WITH EMISSIVITY FOOT**2 
THRMPR(51-54) 0.15 

REFERENCE TEMPERATURE FOR STANDARD LENGTH ZLO DEG R 
THRMPR(81-84) 500.0 

INITIAL TEMPERATURE OF ELEMENT 
THRMPR(91-94) 500.0 

INTEGRATION BOUNDS 
THRMPR(2) 1.0D-8 1.0D-10 


Plotting Locations for Element Temperatures 


KPLOTS(33 1-340) 


Element average temperature 
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ELEMENT DAMPING COEFFICIENT 


Fortran 

Symbol 

lOMKDM(l-lO) 


0MKDMP(3,10) 


Element Twist 


Fortran 

Symbol 

ITWIST 


NKT(IO) 


ZA(10) 


Preset 

Type Value Description 


1*4 0 Control word to invoke use of reference 

frequency 

IOMKDM(K) 0 Reference 

Frequency not used 

IOMKDM(K) 1 Reference 

Frequency used 

R*4 0.0 DO Reference frequency (rad/sec) for 

calculation of model damping coefficient 


(Torsion) Option 


Math 

Symbol Description Units 

Control word to include or exclude twist N.D. 

motion for flexible element. 

ITWIST - 0, twist motion excluded 
ITWIST *= 1, twist motion included 
(preset = 0) 

Set to the number of twist modes desired N. D. 
for a particular flexible element; i. e. , 

NK(K) - 0, 1, 2 or 3 for the k 1 * 1 element. 

Cross-sectional area for the k th flexible in 2 

element, (preset - 3. 757 D-2) 
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Symbol Symbol Description Units 

I2OVI3(10) I 2 /I 3 Ratio of cross-section moments of inertia. 

/ ZjdA 

I 2 /I, (Preset - 1.0) 

/ ZfdA 


ZDQ(10) 


ZJ(10) 


D2(10) 

D3(10) 



f <t>l dA for the cross section. 0 W Is 
warping function, (preset = 7. 463 D-8, 
elliptic cross section) 

in 6 

J 

Torsional constant. 

in 4 


’■/[(‘■•SJ' (-■••!)']■* 



(preset = 2. 118 D-4, elliptic cross 
section) 


d 2 

/Z 2 0 w dA for cross section, (preset = 
0. ODO) 

in 5 


/Z,4> dA for cross section, (preset = 
0. ODO) 

in 5 


CW(3, 10) 
CDW(3, 10) 
CDTW(3, 10) 


Twist angle, (preset = 0. ODO) 

deg 

Twist velocity, (preset = 0. ODO) 

deg/ sec 

Twist damping coefficient. 

ft-lbs 

(preset = 0. ODO) 

rad/sec 


TIP MASS ROTARY INERTIA OPTION 


Fortran 

Math 



Symbol 

Symbol 

Description 

Units 

ITPROT 


Control word to include or exclude tip 
mass rotary inertia. 

ITPROT = 0, rotary inertia excluded 

ITPROT = 1, rotary inertia included 
(preset = 0) 

N.D, 
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Symbol 


Math 

Symbol 


Description 


Units 


NUMTIP(IO) 


Set to 1 for tip mass rotatory inertia 
simulation for the k* element, 
(preset = 0) 

N.D. 

TIPINR(3, 10) 


Principal rotatory inertias of tip mass 
about undeformed element axis. 
TIPINR(I, K) is inertia about I th element 
axis for the k 1 * 1 element tip mass. 

slug ft 

LIBRA TION DAMPER OPTION 


Fortran 

Symbol 

Math 

Symbol 

Description 

Units 

IDAMP 


Control word to exclude or allow 
libration damper motion - 
(Preset value = 0) 

Integer, 

N.D. 



IDAMP = 0 No damper motion 
IDAMP = 1 Damper motion 


1DPHLD 


Control word for simulation of constant 
relative angular velocity of libration dam- 
per mas 0 , or system of booms. 




IDPHLD = 0 No simulation (preset) 




IDPHLD = 1 Constant angular velocity for damper 

NDAMPR 


The number of damper boom elements 
composing the libration damper. 

Integer, 

N.D. 

PHIS 


Stop angle P>r libration damper, 
(preset = 35.0) 

deg 

PHILD 

^LD 

The angular deflection of the libra- 
tion damper boom relative to its 
equilibrium position. Also Euler 
angle in definition of libration damper 
frame (Z) with respect to body Y 
frame, (preset = 0.0) 

deg 

DPHILD 

0LD 

Angular velocity of libration damper 
boom relative to the body, (preset = 0.0) 

deg/ sec 
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Symbol 

Math 

Symbol 

Description 

Units 

BETLD 

*LD 

Euler angle of libration damper Z 
frame with respect to the body 
Y frame, (preset = 0.0) 

deg 

GAMLD 

?LD 

Euler angle of Ubration damper Z 
with respect to body Y frame, 
(preset = 0.0) 

deg 

YIZM(1, n)2 

1 ^emf 

Coordinates of the origin of the 
libration damper Z frames as 

feet 


defined in the body Y frame. 

(preset = 0.0) 

NOTE 

YIZM locates the point of rotation of damper in the main body frame. 

ZBZ locates the root of tne elements with respect to the point where 
rotation takes place (which is specified by the YIZM array) 
for those elements composing the damper. In short, the ZBZ 
array locates the element roots with respect to whatever frame 
you are in. 

For non damper elements use ZBZ array only. 

For damper elements use Z3Z + YIZM 
(one vector for each boom) 

NOTE 

The 2 after the parentheses indicates that it is the second coordinate 
which varies. Also, since ZLO, ZL1, ZLA, and the ZBZ's are i. out 
for both fixed elements and damper booms , the order in which their 
values are punched on an input card is important. The order is: 
data for element booms first, followed by data for damper booms. 

For example, if the user were to set NELMTS 1, NDAMPR 2, 
(implying the spacecraft has a total of three appendages) und the fol- 
lowing card is also input; ZBZ (1,1)2 5.0 7.0 8.0 the program 
would interpret this as the "1" axis of the one rigidly attached ele- 
ment is located 5 feet from the origin of the body frame, the "1" axis 
of the 1st libration damper frame is located 7 L,et from the origin of 
the body frame, and the "1" axis of the 2nd libration damper frame is 
located 8 feet from the origin of the body frame. 
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Symbol 

Description 

Units 

DIN(k, j)l 

Component of the k 111 element (libration damper) 
tip displacement in the j 111 bending mode as 
measured in the element frame (Z) along the 3 
axis of the element frame, (preset = 0.0) 

feet 

DINDOT(k, j)l 

Component of the k 1 * 1 element (libration damper) 
tip velocity in j 1,1 bending mode as measured in 
the element frame (Z) along the 3 axis of the 
frame, (preset = 0.0) 

ft/sec 

DOUT(k, j)l 

Component of the k th element (libration damper) 
tip displacement in the j* bending mode as 
measured in the element frame (Z) along the 2 
axis of the frame . (preset = 0.0) 

feet 

DOUTDT (k, j)l 

Component of the k th element (libration damper) 
tip velocity in the bending mode as measured 

in the element frame (Z) along the 2 axis of the 
frame, (preset = 0.0) 

ft/sec 


NOTE 


In the above arrays, the "1" following the parentheses indicates that 
the first integer in the parenthesis varies as the numbers are read off 
across the card. For example, the card: 

DIN(k, 2)1 1.0 3.0 4.5 2.0 

would be interpreted by the program as: 

DIN(1, 2) = 1.0 
DIN (2, 2) = 3.0 
DIN (3, 2) = 4.5 
DIN(4, 2) - 2.0 etc. 

The arrays DIN, DINDOT, DOUT, and DOUTDT are both input and 
output of the program. Thus if the user wished to start a problem 
with initial tip deflections and velocities for the libration damper 
elements, these quantities would be input. Otherwise they are in- 
ternally set to zero at the ’Start of the problem and are output only. 
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Symbol 

Math 

Symbol 

Description 

Preset 

Value 

Units 

ZK1D 

k lD 

Torsional spring constant for libration 
damper. 


ft-lb/rad 

ZK2D 

k 2D 

Stop spring constant for libration damper. 

2.0 

ft- lb/ rad 

ZMDO 

M ko 

Saturation moment for magnetic 
hysteresis libration damper. 

1.02 D-3 

ft- lb 

ZMDBO 


Initial moment for magnetic hysteresis 
libration damper. 

0.0 

ft-lb 

DECAY 

a 

Exponential decay factor for magnetic 
hysteresis libration damper. 

20.0 

N.D. 

DPRMI(1, 

j)2 

Moments of inertia of the hub of the 
libration damper about the Z frame 
axes of the damper. 

0.0 

slug ft 2 

CNV 

C n v 

Damping coefficient for viscous 
damping by libration damper. 

0.0 

ft-lb/ 
rad/ sec 


EARTH ORIENTED SATE LUTE OPTION - ATTITUDE AND BODY RATES 


Fortran 

Symbol 

Math 

Symbol 

Preset 

Description Value 

Units 

ALFAE 

a 

Roll Angle. The inter- 
mediate Euler angle in a 
2-1-3 rotation of the local 
vertical to body frame. 

0.0 

deg 

BETAE 

& 

Pitch Angle. The first Euler 
angle in a 2-1-3 rotation of the 
local vertical to body frame. 

0.0 

deg 

GAMAE 

7 

Yaw Angle. The last Euler 
angle in a 2-1-3 rotation of the 
local vertical to body frame. 

0.0 

deg 
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Math 


Preset 



Symbol 

Description 

Value 

Units 

OMBC (i) 

(i) = 1. 2, 3 

CJi 

1 B/C 

Components of satellite 
angular velocity with re- 

0.0 

deg/sec 


spect to the local vertical 
frame , with components 
expressed in the body 
frame. 

IBTEST Terminates case (within 0 N.D. 

stack) when BETAE goes 
from - to +. Used in boom 
deployment. (Preset value = 0) 

IBTEST = 0 By-passes option 

IBTEST = 1 Employs option 

BOOM DEPLOYMENT FROM SUN PULSE OPTION 


Fortran 

Symbol 

Description 

Preset 

Value 

ISDPLY 

ISDPLY = 0 By-passes subroutine 
ISDPLY = 1 Employs subroutine 

0 

ISAXIS 

Axis of rotation of spacecraft . (Should be either 1 , 2 , or 3) 

3 

NCROSS 

Number of sun crossings before deployment begins. 

5 

STANG 

Delay angle, (deg) 


ANGTOL 

Angle tolerance, (deg) 


NPRINT 

Number of print-outs from (NCROSS-1) crossings 
to NCROSS (the last) crossing. 

2 

IRAXIS 

Spacecraft body axis upon which sun sensor is 
located. 

3 
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SPINNING BODY OPTION ATTITUDE AND BODY RATES 


Fortran 

Symbol 

Math 

Symbol 

Description 

Units 

PSI1 

* 

First Euler angle in a 3-1-3 rotation 
from local inertial to body frame, 
(preset = 0.0) 

deg 

THETl 

e 

Intermediate Euler angle in a 3-1-3 
rotation from local inertial to body 
frame, (preset = 0.0) 

deg 

PHI1 

* 

Final Euler angle in a 3-1-3 rotation 
from local inertial to body frame, 
(preset = 0.0) 

deg 

OMEG(i) 

(i) = 1, 2, 3 

W iB 

Components of satellite angular ve- 
locity with respect to inertial space, 
with components expressed in body 
frame. Spin axis is the third body 
axis, (preset = 0.0) 

deg/ sec 

ETTA 

ZETTA 

Tf (2nd angle) 
t (1st angle) 

Angles which specify a preferred 
inertial direction of the body spin 
axis with respect to inertial space. 
This is a 3-2 rotation from the 
local inertial frame, (preset = 0.0) 

deg 

deg 

NUTATION VISCOUS RING 

DAMPER OPTION 


Fortran 

Symbol 

Math 

Symbol 

Description 

Units 

IVISCS 


Control word to employ viscous 
ring nutation damper into 
program. (Preset value = 0) 

Integer, 

N.D. 



IVISCS = 1 Employs subroutine 
IVISCS = 0 By-passes subroutine 


VDSCTY 

7 

Kinematic viscosity of liquid in 
nutation damper, (preset = 3.0) 

Centi stokes 

RADTBE 

& 

Radius of nutation damper tube . 

inches 
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Math 

Symbol 

Description 

Units 

RADRNG 

b 

Ra s of nutation damper ring, 
(preset = 10.0) 

inches 

DENSTY 

P 

Density of fluid in nutation 
damper, (preset = 56.16) 

lb/ft 3 

ETTAD 

ZETTAD 

% (2nd angle) 
? D (1st angle) 

Angles which specify the axis of 
the ring damper with respect to 
the body. This is a 3-2 rotation 
from the body frame to the dam- 
per axis. The Y 3 body axis 
initially coincides with the 
damper axis, (preset - 90.0) 

deg 

deg 

YARRAY 

Y,, Y 2 , Y 3 

Components of torque exerted 
by the liquid upon the satellite , 
directed along the damper ring 
axis, (preset = 0.0) 

ft-lb 

OMEGL 

u i 

Initial spatial average angular 
velocity of the liquid relative 
to the ring, (preset = 0.0) 

deg/sec 

ATTITUDE CONTROL MOMENT OPTION 


Fortran 

Symbol 

Math 

Symbol 

Description 

Units 

IATTDE 

Control word to call the attitude control 
subroutines . 

Integer, 

N.D. 


IATTDE - 0 By-pas oes subroutine 
IATTDE = 1 Employs subroutine 
(Preset = 0) 


DTMXA Time increment from initial problem seconds 

time to when the attitude control sys- 
tem is activated (time reference at 
beginning of stack for stacked cases) 

(preset = 1.0) 
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Fortran Math 

Symbol Symbol 


Description 


Units 


PX1 B Fraction of spin period after sun line N.D. 

reference point established when 
control torque is applied to Y, body 
axis. Absolute value should be greater 
than one print interval, (preset = 0.4) 

PXO C Fraction of spin period after sun line N.D. 

reference point established when 
control torque is removed from Y, 
body axis, (preset = 0.5) 

CMX M. Control moment along Y, body axis. ft- lb 

>y Input as a positive number means mo- 

ment applied about positive Y { body 
axis. Input as a negative number 
means moment applied about negative 
Yi body axis, (preset = 0.0) 

NPULSE The number of sequential moment Integer, 

pulses applied to the body once the N.D. 

control system is activated. The 
application rate is one pulse per spin 
period, (preset = 4) 

Constant Torque Levels About Body Axes 


Fortran 

Math 



Symbol 

Symbol 

Description 

Units 

ITORK 


Control word to apply torque about body axes. 
(Preset value = 0) 

ITORK = 0 No torque applied 

ITORK = 1 Torque applied 

Integer, 

N.D. 

CMTORK(3) 


Torque magnitude applied to body axes . 
(preset = 0.0) 

ft— lb 
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SPIN AXIS MOMENT OPTION 


Fortran 

Symbol Description 

ESPEN3 Control word to call spin axis moment subroutine. 

(Preset value = 0) 

ISPIN3 = 0 By-pass subroutine 

1SPIN3 = 1 Employ subroutine 

DTZMA Time increment from initial problem time to 

leading edge of spin moment. 

PZDT Duration of spin moment. 

CMZO Magnitude of moment about the spin axis. Input as 

a positive number means moment applied about the 
positive Y 3 body axis. 

THRUST LOADING OPTION 


Fortran 

Math 


Symbol 

Symbol 

Description 

IPULSE 

N.A. 

Control word to act.. "*e thrust application 

IPULSE = 0 No thrusting 

IPULSE - 1 Apply thrust once 

IPULSE > 1 Apply thrust IPULSE 
times (only if ISPLSE =£ 0. ) 

ISPLSE 

N. A. 

Control word to activate sun crossing time 
to start the thrusting 

ISPLSE = 0 Sun crossing not used 

ISPLSE = 1 Sun crossing used 

ISPNP 

N.A. 

Control word to print out the orbit update 
message (only if ISPLSE = 1, IPULSE >1) 
i.e. , if ISPNP = 5, the orbit update message 
will be printed at every 5th pulse. 

1PLPR1 

N.A. 

Control word for number of thrust pulses 
per spin record 


Units 

Integer , 
N.D. 

seconds 

seconds 
ft— lb 


Units 
N. D. 


N. D. 

N. D. 
N.D. 


56 



Fortran 

Symbol 

TVECTR(3,2) 

TLOCAT(3,2) 

TTIMES(4,2) 


TPARAM(4,2) 


REFANG(2) 


Math 

>ymbol 


< D v> 

{*T> 


•l 

*2 

*3 

u 

A 

B 

C 

D 

a R 


Description Units 

1PLPRP = 1 One pulse 
IPLPRP = 2 Two pulses 

Onl> one or two pulses are allowed. (Preset = 1) 

Unit vector defining the direction of the force N.D. 
applied to the body due to thrusting. This vec- 
tor is defined in the body frame, (preset = 0.0, 

0 . 0 , 1 . 0 ) 

Location in the body frame of the point applica- feet 
tion of the force due to thrusting, (preset = 0.0) 

Times to define thrust variation measured from sec 
the problem starting time (see page A- 16). 

TTIMES (1,1) Start of pulse 

TTIMES (2,1) End of exponential rise 

TTIMES (3,1) End of linear thrust 

TTIMES (4,1) End of pulse 1=1 or 2 

Parameters to define thrust variation 

TPARAM (1,1) Coefficient during exponential rise lb 

TPARAM (2,1) Exponential decay constant during sec -1 
exponential rise 

TPARAM (3,1) Coefficient for linear slope lb/sec 

TPARAM (4,1) Exponential decay constant dur- sec -1 
ing exponential decay 

Angular delay from the Yj axis crossing the deg 
sun line to the initiation of the pulse. 
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JET DAMPING OPTION 


Fortran 

Symbol 

Math 

Symbol 

Preset 

Description Value 

Units 

IJTDMP 


Control word to invoke jet 0 

damping option 

IJTDMP 0 No jet damping 
IJTDMP 1 Jet damping 

Integer, 

N.D. 



NOTE: Jet damping implies 

thrusting and therefore 
IPULSE>0 


TANKCG(3) 


Position vector to center of 0.0 DO 
mass of fuel to be expended 
during thrusting. 

ft 

FUELPP(2) 


Fuel to be expended during 
one thrust pulse 

slugs 

FUELM 


Total mass of fuel. (SLUGS) 

slugs 

RGYFL(3) 


Square of the radius of 
gyration of the fuel mass 
about its own center of mass 

ft 2 

ANGULAR MOMENTUM OPTION 


Fortran 

Symbol 

Math 

Symbol 

Description 

Units 

IHCALC 

N. A. 

Control word to calculate and print the 

N. D. 


angular momentum of the system. 
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Symbol 

Symbol 

Description Units 

IHCALC 


IHCALC = 0 No calculation 

(cont'd) 


IHCALC = 1 Calculate and print 

IHREF 

N.A. 

Control word to compute reference direction N. D. 

for the EPSERR from the system initial 
angula rmom entum . 

IHREF = 0 No calculation 

IHREF - 1 Determining reference direction 
from angular momentum. 


MOMENTUM WHEEL OPTION 


Fortran 

Symbol 

IWHEEL 


Math 

Symbol 

Description 

Units 


Control word to call momentum 
wheel subroutine. 

(Preset value = 0) 

IWHEEL = 0 By-passes subroutine 

IWHEEL = 1 Employs subroutine 

Integer, 

N.D. 


XMOMIN(i) I wj 

(i) = l,2, 3 


VMOM(i) 

(i) = 1, 2, 3 


Inertia of wheel about its spin axis. slug ft 2 

For i = 1, wheel is located on Y, 

body axis , i = 2 wheel is located on 

Y 2 body axis and i = 3, wheel is 

located on the Y 3 body axis. 

(preset = 1. 0, 1. 0, 1. 0) 

Angular velocity of the i wheel with deg/sec 
respect to the body . Input as a posi- 
tive number signifies that the angu- 
lar velocity vector points in the direc- 
tion of the positive body axis. Input as 
a negative number signifies that the 
angular velocity vector points in the 
direction of the negative body axis. 

(preset, = 0.0, 0. 0, 0. 0) 
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Symbol 

Math 

Symhol 

Description 

Units 

DVMOMfl) 

{“w} 

Momentum wheel acceleration 
(negative for deceleration). 

I = 1, 2, 3 (preset = 0. 0, 0. 0, 0. 0] 

deg/sec 2 

1 

VSUR(I) 

w ur 

Upper limit of momentum wheel 
speed. I = 1, 2, 3 (preset = 1. 0D6, 

deg/ sec 
1.0D6, 1.0D6) 

VSDR(I) 

K)Lo 

Lower limit of momentum wheel 
speed. I = 1, 2, 3 (preset = -1. 0D6 

deg/sec 

, -1. 0D6, -1. 0D6) 

MAGNETIC 

MOMENT OPTION 


Fortran 

Symbol 

Math 

Symbol 

Preset 

Description Value 

Units 

IMGMTS 


Control word to exclude or 0 
allow magnetic moments. 
IMGMTS 1 uses constant 
spacecraft dipole defined 
by DPMAG; IMGMTS 2 
provides proportional 
magnetic coil spin up; 

IMGMTS 3 provides con- 
stant level magnetic coil 
spin up. 

Integer, N.D. 

MAGFLD 


Number of harmonics included 
in the representation of the 
earth’s magnetic field (preset = 7). 

Integer, N.D. 

G(l, J) 
H(I, J) 


Coefficients of earth’s magnetic 
field. 


COILS(3) 


Peak coil dipole strength 
along body axes. 

For example, for spin 
up 

COILS(l) -1500.0 
COILS(2) 1500. 0 
for spin down 

COILS(l) 1500.0 
COILS(2) -1500.0 
60 
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Fortran Math Preset 

Symbol Symbol Description Value Units 

Body frame components 0-0 pole-cm 

of the dipole moment of 
the satellite. Positive 
input means that the 
north pole of the dipole 
points in the direction of 
the positive body axis. 

For negative inputs, the 
north pole points in the 
direction of the negative 
body axis. 

NOTE 

An IGRF 1965 Geomagnetic Field model includes 80 spherical 
harmonic coefficients is used (REF6). 

SPACECRAFT ACCELERATION OPTION 

Fortran Math Preset 

Symbol Symbol Description Value Units 

IACOMP Control word to activate 0 N.D. 

the acceleration compu- 
tation for both hub and 
element tip accelerations 
(= 1 ). 

IHUBAC Control word to activate 0 N.D. 

the acceleration compu- 
tation and printout for the 
hub (=1). 

ITIPAC Control word to activate 0 N.D. 

the acceleration compu- 
tation and printout for the 
element tip (=1). 


DPMAG(i) 

(i) = 1, 2, 3 


61 


Accelerometer Location Variables 


Fortran 

Math 

Preset 


Symbol 

Symbol 

Description Value 

Units 

NUMHUB 

N 

Number of accelerometers 6 
and points within the hub 
where acceleration vectors 
are to be computed. 

N.D. 

YHUB(I, J)1 

{^hub} 

Coordinates of the jth 0.0 

accelerometer in the Y 
body frame: 

I = coordinate number 
(i - 1, 2. 3) 

J = accelerometer 
number (j = 1 to N) 

feet 

ALFAEA(J) 

a a 

Second Euler angle de- 0.0 

fining the orientation of 
the No. 1 axis (sensitive 
axis) of the jth accelerom- 
eter to the y body frame 
(rotation about y, ). 

deg 

BETAEA(J) 


First Euler angle defining 0.0 

the orientation of the ac- 
celerometer relative to the 
y body frame (rotation about 

y 2 )> 

deg 

GAMAEA(J) 

y a 

Third Euler angle definirg 0.0 

the orientation of the accel- 
erometer relative to the y 
body frame (rotation about 
y 3 ). (3-a-y or 2-1-3 rotation. 

deg 

ZXI*(k,n) 

Vl.O) 

Shape function evaluated at 1. 0D0 for 
Zj =1.0 for kth element n = 1, 2, 3 

and nth mode, n = 1, 2, 3. 

Normalized cantilever beam 
eigenfunction is used for the 
preset value. 
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Symbol Symbol 


description 


Preset 

Value 


Units 


ZXIP*(k,n) (1. 0) Derivatives of shape func- 1.33D0for 

tion evaluated at z, = 1.0 n - 1 

for kth element and nth 4. 53D0 for 
mode, n = 1, 2, 3. Norm- n = 2 
alized cantilever beam 7. 24D0 for 
eigenfunction is used for n = 3 

preset value. 

ZXIPP*(k,n) XJj(l.O) Second derivative of shape 0. 0D0 for 

function evaluated at 1 X = n = 1. 2, 3 
1.0 f--~ kth element and 
nth mode, n = 1, 2, 3. 

Normalized cantilever eigen- 
function is used for the re- 
set value. 

ZZNP*(k,I,n) f z np ] Axial displacement shape 

function evaluated at 
z, =1.0 for kth element. 

I = 1 to n; J = i to n. 

ZZNP should always be 
input as a square matrix, 

F / ST FOURIER TRANSFORM (FFT) ANALYSIS 

Fortran Math 

Symbol Symbol Description Units 

ICSD Control word to activate the FFT analyses 

subroutine. 

ICSD = 0 No FFT analysis 
ICSD «= N FFT analysis activated 
(preset = 0) 

The integer N requests that N frequencies be 
extracted from the selected KPLOTS data sets. 

Also, the following must be input for 
FFT analysis: KPLOTS(253) 

set to 2 for plot and FFT analysis 
set to 3 for FFT analysis only. 

*When other than a normalized c" .itilevel beam shape function is desired, the user should input the appro- 
priate values of ZXI , ZXIP, ZXIPP and ZZNP. 
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NOTE 


The value of N should not exceed the number of frequencies that can teasonably 
be expected to exist in the data. The range of permissible values for N are from 
1 to 10. 1 < N < 10. 

There are two precautions the user should take when employing the Fast Fourier 
Transform analysis option. First, the data points to be analyzed in the Fast 
Fourier Transform are controlled by the values of TSTOP and FREQ, i. e. , 

N = Sample size = Y req 1 

In other words, the data set input to the FFT subroutines is exactly the same as 
that appearing in the normal printout. Therefore, the input TSTOP value should 
at least cover more than the two longest vibrational periods of the data set to be 
analyzed. The value of FREQ should not be too coarse in order that the higher 
harmonics can be picked up by the analysis. Furthermore, it is preferred that 
the data set to be analyzed is in steady state condition. Therefore, it is recom- 
mended that the user first run the program without using the FFT analysis option 
to obtain a general idea about the periodic behavior of the data sets to be anal- 
yzed, and then run the program with the FFT option inputting the appropriate 
TSTOP and FREQ values. Second, because the data sets gener ited bv FSD pro- 
gram are all deterministic, the autocorrelation function analysis may not be nec- 
essary; hence, IAFM(4) should be set equal to 0. 

ATMOSPHERIC DENSITY MODEL OPTION 


Fortran Math 

Symbol Symbol Description Units 


IDRAG Control word to include or exclude atmos- 

phere models. Set to 1 for 1962 standard 
atmosphere. Set to 2, 3 or 4 for Jacchia- 
Roberts atmospheric model (see below). 
Preset = 1 

IJACCT Modified Julian date of first data on 

Jacchia- Roberts file (10300 corresponds 
to March 19, 1969). Note: inp” + value 

is used o- ■ for IDRAG3. 
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Symbol 


Math 

Symbol 


Description 


Units 


KP(21, 8) Three hour magnetic activity indices for 

21 day period. Note: input values are 
used only for IDRAG 2 or 3. 

TC(20) Exospheric temperatures. Note: input 

values used only for IDRAG 2 or 3. 

ER(3) Earth's rotation rate for rotating atmos- 

phere calculation. 

(preset - 0. 0D0 
0. ODO 

7. 2921 ID- 5 rad/sec) 

NOTE 

For IDRAG 2 the modified Julian date for the Jacchia-Roberts 
data file is internally set to the simulation date defined by 
IDATE. For IDRAG 3 the input modified Julian date defined 
by IJACCT is used. For IDRAG 4 the Jacchia-Roberts data is 
to be read from unit 15 (360/95). 

Linear Varying Drag Loading (for tethered mass) 


XLTEST Linear varying drag loading is calculated 

for any element whose length is greater 
than XLTEST. (preset = 1. 0D + 06) 

NOTE 

The arrays for pressure calculations can be input but are 
preset to appropriate values in BLOCK data. 

GROUND SIMULATION OPTION 
Fortran Math 

Symbol Symbol Description 

IGRUND Control word to activate the ground 

simulation environment. 

IGRUND = 1 Activates option 
IGRUND = 0 Does not activate option 
(preset = 0) 


feet 


Units 

NA 
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Symbol 

IGASBR 


IALTUD 


ALFAEG 

DELTAG 

ALTUDE 

OMGY(l-3) 

GACC(l-3) 


Math 

Symbol Description 

Control word to activate the air bearing 
ground simulation. If IGRUND = 1 and 
IGASBR = 0, the spin platform (single 
axis) capability is used. If IGRUND = 1 
and IGASBR = 1, air bearing (3 axis) 
capability is used, (preset = 0) 

Control word to activate the computation 
of the local verticle direction of the 
gravitational acceleration at the ground 
testing site altitude. If IALTUD = 0, 
gravitational acceleration input is 
necessary; if IALTUD = 1, the altitude 
input is necessary. 

Oq First rotation angle of 3-1 rotation from 

local geographical frame to the body frame 
(or the right ascention angle), see page A-19. 
(Preset = 0. 0D0) 

5 g Second rotation angle of a 3-1 rotation 

from local geographical frame to the body 
frame (or the declination angle of spin or 
number 3 body axes in local geographical 
frame). (Preset = 0. 0D0) 


h c 


i w cf 

JAc! 


The altitude of the ground test site. Used 
to compute the gravitational acceleration 
at the altitude when IALTUD = 1. 

(Preset = 0. 0D0) 

Angular velocity vector of body axes 
expressed in the local geographical frame. 
(Preset = 0. 0D0) 

Three components of acceleration vector 
input along the three axes of the local 
geographical frame. A G1 is along the 
east local horizontal, A G2 is along the 
north local horizontal and A G3 is along 
the local vertical direction. 

(Preset = 0. 0D0 
0. ODO 

32. 145(552D0) 
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Units 

NA 


NA 


deg 

deg 

km 

deg/sec 

ft/sec 2 



SECONDARY BODY SIMULATION 


The effect of a rigid secondary body connected to the primary body through 
a three degree of freedom rotary (universal) joint can be simulated. The con- 
necting joint has either one, two, or three degree of freedom selected by in- 
put. The suspension simulation for the secondary body includes a bilinear 
spring, viscous damping, and angular stops (limits) for each axis. 

A ras ering capability of the secondary body is included so that the azimuth 
angle can be swept through at a specified rate: and at the angular stop an 
incremental step change in elevation may be made. The secondary body then 
moves in the opposite azimuth direction to the opposite stop. etc. Rastering 
stops when the number of input elevation steps have been completed. 

Fortran Preset 

Symbol Type Value Description 


I2BDY 

1*4 

0 

Control word to add secondary body 
to simulation 




I2BDY 0 No secondary body 
I2BDY 1 Secondary body present 

NDOF2 

1*4 

0 

Number of degrees of freedom of 
secondary body 




NDOF2 1 Rotation about the 2 

axis of secondary body. 
NDOF2 2 Rotation about the 3 & 2 
axes of secondary body. 
NDOF2 3 Rotation about the 3. 1 & 
2 axes of 5 -rndary body. 

SECM 

R*8 

0.0D0 

Mass of secondary body. v «). 
NOTE. The mass of the secondary 
body must be included in 
the system mass ZMS. 

SECK3.3) 

P*8 

0.0D0 

Mass moments of inertia of secondary- 
body about its own center of mass. 
(slugs-ft 2 ). 

ZBAR2(3) 

R*8 

O.ODO 

Location of the secondary body center 
cf mass measured from the pivot point 
which is the secondary body fixed 
reference frame origin (ft.). 
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SECONDARY BODY SIMULATIONS (Cont.) 


Fortran 

Symbol 

Type 

Preset 

Value 

Description 

YI02(3) 

R*8 

0.0D0 

Location of the pivot point for the sec- 
ondary body measured in the main body 
reference frame, (ft.). 

SBKK3) 

R*8 

O.ODO 

Spring constant for secondary body sus- 
pension. Applies to angles less than the 
stop angle, (ft-lb/rad). 

SBK2(3) 

R*8 

0.0D0 

Spring constant for secondary body sus- 
pension. Applies to angles greater than 
the stop angle, (ft-lb/rad). 

SBSTA(3) 

R*8 

O.ODO 

Stop angle for secondary body suspen- 
sion. (rad). 

SBDM(3) 

R*8 

O.ODO 

Damping coefficients for secondary body 
motions, (ft-lb-sec/rad). 

GAM 20 

R*8 

O.ODO 

Rotation about the three axis in a 3-1-2 
angle set from the main body frame to the 
secondary body equilibrium frame, (deg). 
NOTE: The secondary body equilibrium 
frame is that orientation at which 
spring have zero restoring torque. 

ALP20 

R*8 

O.ODO 

Rotation about the one axis in a 3-1-2 
angle set from the main body frame to 
the secondary body equilibrium frame, 
(deg). 

BET20 

R*8 

O.ODO 

Rotation about the two axis in a 3-1-2 
angle set from the main body frame to 
the secondary body equilibrium frame, 
(deg). 

GAM21 

R*8 

O.ODO 

Initial condition for rotation about the 
three axis of the secondary body. Mea- 
sured from the equilibrium frame to the 
secondary body fixed frame, (deg). 

ALP21 

R*8 

O.ODO 

Initial condition for rotation about the 


one axis of the secondary body. Mea- 
sured from the equilibrium frame to the 
secondary body fixed frame, (deg). 
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SECONDARY BODY SIMULATION (Cont.) 


Description 


Fortran Preset 

Symbol Type Value 

BET21 R*8 0.0D0 


OM2K3) 

R*8 

0.0D0 

SBUP(2) 

R*8 

2*1.0D-3 

SBDN(2) 

R*8 

2*I.OD-5 


Initial condition for rotation about the 
two axis of the secondary body. Mea- 
sured from the equilibrium frame to the 
secondary body fixed frame, (deg). 

NOTE: The motion of the secondary 

body is simulated in terms of a 
3-1-2 set of angles from the equi- 
librium frame to a secondary body 
fixed frame. If only two degrees 
of freedom are specified, the sim- 
ulation uses a 3-2 set of angles. 

If one degree of freedom is speci- 
fied the simulation uses a rotation 
about the 2 axis. 

Initial conditions for angular rates of the 
secondary body about the secondary body 
fixed axes, (deg/sec). 

Integration upper bounds for secondary 
body angles and angular rates. 

Integration lower bounds for secondary 
body angles and angular rates. 


Plotting locations for secondary body variables 


Fortran 

Symbol 

KPLOTS (236) 
KPLOTS (237) 
KPLOTS (238) 


Description 

Rotation of secondary body about 3 axis 
t sb 

Rotation of secondary body about the carried 1 axis 
a SB 

Rotation of secondary body about the carried 2 axis 
0SB 
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Plotting locations for secondary body variables (Cont.) 


Fortran 

Symbol 

KPLOTS (239) 
KPLOTS (240) 

KPLOTS (241) 

KPLOTS (242) 


KPLOTS (243) 


KPLOTS (244) 


Description 


Relative angular rate of secondary body about 3 axis 

Relative angular rates of secondary body about the carried 1 
axis of secondary body 

*SB 


Relative angular rates of secondary body about the carried 2 
axis 

£sb 

Component of secondary body relative angular velocity on 1 
axis of secondary body 

W 1SB 

Component of secondary body relative angular velocity on 2 
axis of secondary body 

W 2SB 

Component of secondary body relative angular velocity on 3 
axis of secondary body 

W 3SB 


Secondary Body Rastering Input 


Fortran 

Symbol 

Type 

Preset 

Value 

IRAST 

1*4 

0 

IARST(3) 

1*4 

0 


Description 

Control word to invoke prescribed rastering motions 
for secondary body. 

IRAST = 0 No rastering 

IRAST = 1 Rastering prescribed 

Control word to specify type of rastering cycle to be 
invoked on each axis. 


IARST(1) = 0 No motion 
IARST(I) = 1 Motion of Type 1 
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Preset 

Type Value Description 

IARST(I) = 2 Motion of Type 2 
1=1 to J 

I = 1 motion about the 3 axis 

1 = 2 motion about the 1 axis 

I = 3 motion about the 2 axis 


NOTE: The types of motion cycles are shown in 
Figure 20 and 21 on pages A-20, 21. 


1RSCY(3) 

1*4 

0 

Number of cycles of rastering motion to be per- 
formed on each axis. 

DELA(3) 

R*8 

0.0D0 

Angular motion to be performed during a single 
cycle of rastering motion for each axis. (deg). 

TAUA(4,3) 

R*8 

0.0D0 

Time increments to define rastering cycle for 
each axis. See Figure 20 and 21. 


TAUA (1,1) = T, 

TAUA (2, 1) = T 2 
TAUA (3, 1) = T 3 
TAUA (4, I) = T 4 

ANG20(3) R*8 0.0D0 Initial angle for the start of the rastering cycle. 

(deg). 

ACTUATOR INITIATION FROM ZERO CROSSING OF FILTERED STATE VARIABLES 

Thrusters and momentum wheels can be activated from either positive or negative 
going zero crossings of spacecraft state variables such as accelerometer readings, 
body components of spin vector, body components of earth’s magnetic field, and 
body components of ihe sun vector. A fixed time delay is permitted by input 
starting from the time zero crossing of the state variable. Also, a simple filter of 
the form: 

Kts 

(1 + T S) 2 


Fortran 

Symbol 
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where K = gain factor, input 
s = Laplace operator 
r ~ time constant, input 

is included to take out any d.c. bias thai may be present in the state variable. 
The momentum wheel speed is driven in a triangular wave form between limits 
(VSUK and VSDR) at the frequency specified by input. 


Fortran 

Symbol Type 

Preset 

Value 

Descripiion 

IACFLT(20) R*8 

0 

Control integers for filtered sensor signal to 
actuate momentum wheel cycling or pulsed 
thrusting. 

IACFLS(l) = 0 
lACFLT(l) = 1 


No filtered sensor simulation. 
Filtered sensor simulation. 

IACFLTC) = 1 
IACFLTC) = 2 


Pulsed thrusting to be simulated. 
Momentum wheel cycling to be simulated. 

1ACFLTO) = 1 
IACFLK3) = 2 
IACFLK3) = 3 
IACFLK3) = 4 


Sensor measures magnetic field. 

Sensor measures an acceleration component. 
Sensor measures body angular velocity. 
Sensor measures solar intensity. 

IACFLK4) = 0 


Use natural initial conditions for filter inte- 
grates. 

IACFLT(4) = 1 


Use input initial conditions for filter inte- 
grater. 

ACPARM(20) R*8 

0.0D0 

Numerical data required for filtered sensor 
signal simulation. 

ACPARM(l) 


Filter gain. 

ACPARM12) 


Filter time constant. 

ACPARM(3-5) 


Unit vector defining sensor direction in the 
body frame. 

ACPARM(6-8) 


Sensor position in the body frame, (for ac- 
celeration sensor). 

ACPARM(9) 


Cycle half period for momentum wheel 
cycling. 
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Actuator Initiation from Zero Crossing of Filtered State Variables (Cont.) 

Fortran 

Symbol Type 

Preset 

Value 

Description 

ACPARM(IO) 


Time delay for initiation of actuator after 
zero crossing. 

ACPARM(17) 


Initial condition for filter integrates 

ACPARM(18) 


Initial condition for filter integrates 

ACPARM(19) 


Upper bound for filter irtegrater. 

ACPARM(20) 


I ower bound for filter integrates 

Plotting locations for filter state variables 

Fortran 

Symbol 


Description 

KPLOTS1245) 


Sensor output 

KPLOTS(246) 


Filter output 

KPLOTS(247) 


Filter first integrater output 

DUAL-SPIN SPACECRAFT CONTROL SYSTEM (DE-B) 


The input Fortran symbols for the DE-B control and active nutation damper 
system are given below. A block diagram of the control system is given in 
Appendix A, page A-22 along with some mathematical description of the for- 


mulation. 




Fortran 

Symbol 

Type 

Preset 

Value 

Description 

KNTRL(IO) 

1*4 

0 

Vector of control integers for DE-B control 


system simulation 


KNTRL(l) = 0 No control system 
KNTRL(l) = 1 Control system with second 
order sensors 

KNTRL(l) = 2 Control system with fourth 
order sensors 

(2) KNTRL(2) = 0 No nutation damper 

KNTRL(2) = 1 Primary damper circuit 


KNTRL 

(1) 
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Fortran 

Symbol 

Type 

Preset 

Value 

Description 

(2) 



KNTRL(2) = 2 Offset pointing damper circuit 

(3) 



KNTRL(3) = 0 No 

KNTRL(3) = 1 First order tachometer 

(4) - (8) 



Not used 

(9) 



KNTRL(9) = 1, I is starting integer for noise 
generator. I must be odd and should have 6 
or 7 digits. 

(10) 



KNTRL(IO) = 0 No noise channel? 

KNTRL(IO) = 3 Noise generated for sensors & 
bial voltage. Use only 0 or 3. 

CPARM(43) 

R*8 

(1-30) 0.000 
(31-40)1.000 
(41-43)0.000 

Control system parameters for DE-B control 
system simulation 

CPARM 

Math 

Symbol 

Units 

Description 

(1) 

T s 

sec. 

Sensor time constant 

(2) 

T ! 

sec. 

Lead term in pitch compensation 

(3) 

r 2 

sec. 

Lag term in pitch compensation 

(4) 

t F 

sec. 

Tachometer time constant 

(5) 

K s 

volts/rad. 

sensor gain 

(6) 

K c 

volts/volt 

Pitch amplifier gain 

(7) 

K a 

volts/volt 

Power amplifier gain 

(8) 

Kf 

volts/(rad/sec) Tachometer gain 

(9) 


volts/(rad/sec) Motor back EMF constant 

(10) 

K « 

ft.lbs/volt 

Motor torque constant 

(11) 



Not used for input 
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CPARM 

Math 

Symbol 

Units 

Description 

(12) 

V lin, 

volts 

Voltage limit in compensation amplifier 

(13) 

Si 

volts/volt 

Gain in primary damper circuit 

(14) 

Vi 

sec. 

Time constant in primary damper circuit 

(15) 

v b 

volts 

Bias voltage 

(16) 

r co 

ft-lbs 

Coulomb friction torque 

(17) 

ft 

min 

rad /sec 

S 2 

co 

Test relative wheel speed to avoid coulomb 
friction torque discontinuity at zero speed 


T c ~~ T c0 

fynin + l fi cj 


(18) 

Ps 

volts/voits 

In fourth order sensor. See p. A-22. 

(19) 

S2 

sec. 

Gain in offset pointing damper circuit 

(20) 

7 )x2 

N.D. 

Time constant in offset pointing damper 
circuit 

(21) 


N.D. 

Sign of dampei circuit output set 1.0 
or -1.0 

(22) 


volts 

Roll sensor output limit 

(23) 


ft. -lbs. 

Motor torque output upper limit 

(24) 


ft. -lbs. 

Motor torque output lower limit 

(25-30) 



Not used 

(31-33) 



Noise model SIGMA for pitch, roll and 
vo. age bias respectively 

(34-35) 



Not used 

(36-38) 



Noise model LAG for pitch, roll and vol- 
tage bias respectively 

(39-40) 



Not used 


75 



CPARM 

Math 

Symbol 

Units 

Description 

(41) 


rad. 

Pitch sensor bias 

(42 \ 


rad. 

Roll sensor bias 

(43) 



Not used 

Fortran 

Symbol 

Type 

Preset 

Value 

Description 

SVCS(20) 

R*8 

0.000 

Initial conditions for control system state 
vector 

SVCS 


Units 

Description 

(1) 


volts 

Pitch sensor output 

(2-4) 


volts 

Pitch sensor dynamics 

<5) 



Not used 

(6) 


volts 

Roll sensor output 

(7-9) 


volts 

Roll sensor dynamics 

(10) 



Not used 

(11) 


volts 

State variable for pitch compensation 
amplifier 

(12) 



Not used 

(13) 


volts 

Tachometer output 

(14) 



Not used 

(15) 


rad/sec 

Wheel speed 

(16-18) 



Not used 

09) 


volts 

Nutation damper 

(20) 


volts 

Nutation damper 
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Fortran 


Preset 


Symbol 

Type 

Value 

Description 

GN1C(10) 

R*8 

0.0D0 

Initial conditions for noise model channels 


GNIC(l) Pitch channel 
GNIC(2) Roll channel 
GN1C(3) Bias voltage channel 
GNIC(4-10) Not used 

CSUP(20) R*8 1.013-2 Upper bound on difference bet veen i 'dieted 

and corrected control system state vector. Lo- 
cation in CSUP corresponds to the location of 
the variable in the state vector initial condition 
array SVCS 

CSDN(20) P : 8 1.0D-4 Lower bound cn difference between predicted 

and corrected control system state vector. Lo- 
c?tion in CSDN corresponds to the location of 
the variable in the state vector initial condition 
array SVCS 

In addition to the input described above, the following m> st be given: 

1WHEEL 1 Makes call to momentum wheel subroutine (WHEELS) 

XMOMIN(2) — Inertia of momentum wheel (slug— ft 2 ) 

It is recommended that integration control bounds oe given for critical control 
variables to prevent integration time step from exceeding the time constants of 
the closed-loop system. For example, if fastest component in the control system 
has a time constant of .1 sec., it is unreasonable to expect accurate simulation re- 
suits with a larger time step. Setting the integration bounds to some reasonable 
(small) fraction of the nominal value (eg 10“ 2 to 10‘ 4 ) will ensure that the inte- 
gration errors are of the same order (e.g. 2 to t p’ace accuracy). The time stop 
corresponding to this level of accuracy will be determined internally. 

The output for the DE-B control system simulation includes both printed data 
and plots. The printed output is as follows: 


Fortran Symbo 1 

Description 

Units 

PTCH Out 

Pitch channel sensor output 

volts 

ROLL Out 

Roll channel sensor output 

volts 
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Fortran Symbol 

Description 

Units 

COMP Out 

Output of compensation amplifier 

volts 

TACH Out 

Output of tonometer 

volts 

TMOTOR 

Torque output of momentum wheel motor 

ft-lbs 

WHL SPD 

Momentum wheel speed 

rad/sec 

NUTD Out 

Nutation damper phase shift circuit output 

volts 

The output available for plotting is the entire state vector for the 
The locations and definitions of these variables are as follows. 

control system. 

K PLOTS 
Address 

De— : "tion 

Units 

216 

Pitch senior output 

volts 

217-219 

Pitch sensor dynamics 

volts/sec, etc. 

220 

N 't used 


221 

Roll sensor output 

volts 

222-224 

Roll sensor dynamics 

volts/sec, etc. 

225 

Not used 


226 

Output of pitch compensation amplifier 

volts 

227 

Not used 


228 

Tachometer output 

volts 


TWO AXIS GIMBLE SIMULATION (No Active Control) 


Fortran 

Symbol 

Type 

I reset 
Value 

Description 

IGMBL 

1*4 

0 

Control word for two axis gimble simulation 
IGMBL 0 No gimble simulated 
IGMBL 1 Gimble simulated 

AZIN(3,3) 

R*8 

0.0D0 

Moments of inertia of the azimuth platform about 
its own center of mass (slug-ft 2 ) 

AZAX(3) 

R*8 

0.0D0 

Position vector in the body frame to the origin of 
the gimble motion reference frame. This position 
is a point on the azimuth motion axis, (ft) 

AZCG(3) 

R*8 

O.ODO 

Position vector to the azimuth platform center of mass 
in the azimuth motion reference frame, (ft.) 

AZMS 

R*8 

0.0D0 

Mass of the azimuth platform (slugs) 

ELIN(3,3) 

R*8 

O.ODO 

Moments of inertia of the elevation platform about 
its own center of mass, (slug-ft 2 ) 

ELAX(3) 

R*8 

O.ODO 

Position vector in the azimuth motion reference 
frame. This position is a point on the elevation 
motion axis, (ft) 

ELCG(3) 

R*8 

O.ODO 

Position vector to the elevation platform center 
of mass in the elevation motion reference frame, (ft) 

ELMS 

R*8 

O.ODO 

Mass of the elevation platform (slugs) 

GMK1(2) 

R*8 

O.ODO 

Spring constant for two axis gimble suspension. 
Applies to angles less than the stop angle, (ft-lb/rad) 




GMK1(1) Azimuth Axis 
GMK1(2) Elevation Axis 

GMK2(2) 

R*8 

O.CDO 

Spring constant for two axis gimble suspension. Applies 


to angles greater than the stop angle, (ft-lb/rad) 
GMK2(1) Azimuth Axis 
GMK2(2) Elevation Axis 


GMDMP(2) R*8 0.0D0 Viscous damping constant for two axis gimble 

suspension, (ft-lb-sec/rad) 

GMDMP(1) Azimuth Axis 
GMDMP(2) Elevation Axis 
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TWO AXIS GIMBLE SIMULATION (No Active Control) CONT. 


Fortran 


Preset 


Symbol 

Type 

Value 

Description 

GMSTP(2) 

R*8 

0.0D0 

Stop angle for two axis gimble suspension, (rad) 
GMSTP(1) Azimuth Axis 
GMSTP(2) Elevation Axis 

AZIMO 

R*8 

0.01)0 

Rotation about the three axis in a 3-1-2 angle 
set from the main body frame to the gimble motio 
reference frame, (deg) 

NOTE: The gimble motion reference frame is a 
main body fixed ' rame 

ROLLO 

R*8 

0.0D0 

Rotation about the one axis in a 3-2-1 angle set 
from the main body frame to the gimble motion 
reference frame (deg) 

ELEVO 

R*8 

0.0D0 

Rotation about the two axis in a 3-1-2 angle set 
from the main body frame to the gimble motion 
reference frame (deg) 

AZIMI 

R*8 

O.ODO 

Initial azimuth angular position for the azimuth 
platform (deg) 

AZIMID 

R*8 

O.ODO 

Initial azimuth angular velocity for the azimuth 
platform (deg/sec) 

ELEVI 

R*8 

O.ODO 

Initial elevation angular position for the elevation 
platform (deg) 

ELEVID 

R*8 

O.ODO 

Initial elevation angular velocity for the elevation 
platform (deg/sec) 

GMUP(2) 

R*8 

1.0D-3 

Integration upper bounds for gimble angles and 
angular rates 

GMUP(1) (rad) 

GMUP(2) (rad/sec) 

GMDN(2) 

R*8 

1.0D-5 

Integration lower bounds for gimble angles and 
angular rates 

GMDN(l) (rad) 

GMDN(2) (rad/sec) 
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TWO AXIS GIMBLE SIMULATION (No Active Control) CONT. 


Plotting locations for two axis gimble motion variables 


Fortran 

Symbol Description 

KPLOTS(297) Azimuth angle for the gimble azimuth platform 

KFLOTS(298) Azimuth angular rate for the gimble azimuth 

platform 


KPLOTS{299) 


Elevation angle for the gimble elevation platform 


KPLOTS(300) 


Elevation angular rate for the gimble elevation 
platform 


(LEF BLANK) 
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CONTROL SYSTEM FOR TWO AXIS PLATFORM MAGNETIC TRACKING 


Fortran 

Symbol Type 
IPLTCSO-20) 1*4 

IPLTCS(l) 


IPLTCS(2) 


IPLTC3(3) 

IPLTCS(4) 

IPl.TCS(5-20) 
PCSPRM(IOO) R*8 
PCSPRMO-9) 
PCSPRM(IO) 
PCSPRM(ll) 

PCSPRM( 1 2) 


Preset 

Value Description 

0 Control words for gimble platform control system 

simulation 

Control word to invoke the two axis platform control 
system option 

1PLTCS(1) 0 No control system 
IPLTCS(l) 1 Control system 

NOTE : If a control system is invoked the gimble 
simulation must also be invoked 

Control word to specify magnetometer is mounted at 
the tip of an element 

IPLTCS(2) 0 Magnetometer in body 
IPLTCSt2) K Magnetometer on Kth element 

NOTE. K is defined in terms of the internal 
program element ordering system 

Control word for guassian noise transfer function for 
the three axes of the magnetometer 

IPLTCS(3) 0 No transfer function 
IPLTCS(3) 1 Transfer function 

Number of magnetometer samples to be averaged 
for smoothing of magnetometer output 

IPLTCS(4) N average N samples 
N > 1 

Not used 

0.0D0 Platform control system parameters 
Not used 

Sampling time for PID digital controller (sec) 

Azimuth axis quantatisation level for platform 
position (rad) 

Azimuth PID integrator upper saturation limit (rad) 

Azimuth PID integrator lower saturation limit (rad) 
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CONTROL SYSTEM FOR TWO AXIS PLATFORM MAGNETIC TRACKING (CONT) 


Fortran 

Preset 


Symbol Type 

Value 

Description 

PCSPRM(14) 


Azimuth PID proportional train K p 

PCSPRM(15) 


Azimuth PID integrator gain Kj 

PCSPRM(16) 


Azimuth PID derivative gain Kp 

PCSPRM( 17-20) 


Not used 

PCSPRM(21) R*8 

0.0D0 

Elevation axis quantatisation level for platform 
position (rad) 

PCSPRM(22) 


Elevation PID integrator upper saturation limit 

PCSPRM(23) 


Elevation PID integrator lower saturation limit 

PCSPRM(24) 


Elevation PID proportional gain K p 

PCSPRM(25) 


Elevation PID integrator gain Kj 

PCSPRM(26) 


Elevation PID derivative gain Kp 

PCSPRM( 27-29) 


Not used 

PCSPRM(30) 


Integration upper bound for magnetometer first 
order lag transfer function 

PCSPRM(31) 


Integration lower bound for magnetometer first 
order lag transfer function 

PCSPRM(32) 


Magnetometer one axis bandwidth 

PCSPRM(33) R*8 

0.0D0 

Magnetometer two axis bandwidth 

PCSPRM(34) 


Magnetometer three axis bandwidth 

PCSPRM(35) 


Magnetometer sampling rate (sec) 

PCSPRM(36) 


Computational delay (sec) 

PCSPRM(37-40) 


Not used 

PCSPRM(4 1 ) 


Azimuth platform amplifier gain K A 
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CONTROL SYSTEM FOR TWO AXIS PLATFORM MAGNETIC TRACKING (CONT) 


Fortran 

Preset 


Symbol Type 

Value 

Description 

PCSPRM(42) R*8 

0.0D0 

Azimuth platform motor torque constant Kj 

PCSPRM(43) 


Azimuth platform back EMF constant Kg 

PCSPRM(44) 


Azimuth platform motor torque upper limit Lj 
(ft-lbs) 

PCSPRM(45) 


Azimuth platform motor torque lower limit L-> 
(ft-lbs) 

PCSPRM(46) 


Azimuth platform coulomb friction torque 
constant (ft-lbs) 

PCSPRM(47) 


Azimuth platform minimum angular rate for coulomb 
friction torque (rad/sec) 

PCSPRM(48-50) 


Not used 

PCSPRM(Sl) 


Elevation platform amplifier gain 

PCSPR M(52) 


Elevation platform motor torque constant Kj 

PCSPRM(53) 


Elevation platform back EMF constant Kg 

PCSPRM(54) 


Elevation platform motor .orque upper limit Lj 
(ft-lbs) 

PCSPRM(55) 


Elevation platform motor torque lowe r limit L? 
(ft-lbs) 

PCSPRM(56) 


Elevation platform coulomb friction torque constant 
(ft-lbs) 

PCSPRM(57) R*8 

0.0D0 

Elevation platform minimum angular rate for 
coulomb friction torque (rad/sec) 

PCSPRM(58-79) 


Not used 

PCSPRM(80) 


Amplitude of sinusoidal noise added to one axis 
magnetometer measurement (gauss) 

Amplitude of sinusoidal noise for two axis (gauss) 
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CONTROL SYSTEM FOR TWO AXIS PLATFORM MAGNETIC TRACKING (CONT) 


Fortran 

Symbol Type 

Preset 

Value 

Description 

PCSPRM(82) R*8 

0.0D0 

Amplitude of sinusoidal noise for three axis (gauss) 

PCSPRM(83) 


Phase of sinusoidal noise added to one axis magneto- 
meter measurement (deg) 

PCSPRM(84) 


Phase of sinusoidal noise for two axis (deg) 

PCSPRM(85) 


Phase of sinusoidal noise for three axis (deg) 

PCSPRM(86) 


Frequency of sinusoidal noise added to one axis 
magnetometer measurement (cps) 

PCSPRM(87) 


Frequency of sinusoidal noise for two axis (cps) 

PCSPRM(88) 


Frequency of sinusoidal noise for three axis (cps) 

i'JSPRM(89) 


Not used 

PCSPRM(90-95) 


Used internally. Not input 

PCSPRM(96-1 00) 


Not used 

Plotting locations for platform control system variables 

Input 

Symbol 


Description 

KPLOTS(301) 


Output of first order lag transfer function for magnet- 
ometer one axis 

KPLOTS(302) 


Output of first order lag transfer function for magnet- 
ometer two axis 

KPLOTS(303) 


Outi it of first order lag transfer function for magnet- 
ometer three axis 

KPLOTS(304) 


Azimuth error output 

KPLOTS(305) 


Elevation error output 

KPLOTS(306) 


Azimuth PID digital controller output 
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CONTROL SYSTEM FOR TWO AXIS PLATFORM MAGNETIC TRACKING (CONT) 


Input 

Symbol Description 


KPLOTS(307) 

KPLOTS(308) 

KPLOTS(309) 


Elevation PID digital controller output 
Azimuth platform drive motor torque 
Elevation platform drive motor torque 


TWO AXIS DAMPER GIMBLE SIMULATION (No Active Control) 


Fortran 


Preset 


Symbol 

Type 

Value 

Description 

IGMBLD 

1*4 

0 

Control word for two axis damper gimble simulation 
IGMBLD 0 No gimble simulated 
IGMBLD 1 Gimble simulated 

DELIN (3,3) 

R*8 

0.0D0 

Moments of inertia of the elevation platform about 
its own center of mass (slug-ft 2 ) 

DELAX(3) 

R*8 

J.ODO 

Position vector in the azimuth motion reference 
frame. This position is a point on the elevation 
motion axis, (ft) 

DELCG(3) 

R*8 

0.0D0 

Position vector to the elevation platform center of 
mass in the elevation motion reference frame (ft) 

DELMS 

R*8 

0.0D0 

Mass of the elevation platform (slugs). 

DMK1(2) 

R*8 

0.0D0 

Spring constant for two axis gimble suspension. 
Applies to angles less ihan the stop angle (ft-lb/rad) 

DMK1(2) Elevation axis 

DMK2(2) 

R*8 

0.0D0 

Spring constant for two axis gimble suspension. Applies 
to angles greater than the stop angle, (ft-lb/rad) 
DMK2(2) Elevation Axis 

DMDMP(2, 

R*8 

0.0D0 

Viscous damping constant for two axis gimble 
suspension (ft-lb-sec/rad) 

GMDMP(2) Elevation axis 

DMSTP(2) 

R*8 

0.0D0 

Stop angle for two axis gimble surpension (rad) 
GMSTP(2) Elevation axis 
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TWO AXIS DAMPER GIMBLE SIMULATION (No Active Control) CONT 


Fortran 

Symbol 

Type 

Preset 

Value 

Description 

DLEVI 

R*8 

0.0D0 

Initial elevation angular position for the elevation 
platform (deg) 

DLEVID 

R*8 

0.0D0 

Initial elevation angular velocity for the elevation 
platform (deg/sec) 

DMUP(2) 

R*8 

1.0D-3 

Integration upper bounds for gimble angles and 
angular rates 

DMUP(1) (rad) 

DMUP(2) (rad/sec) 

DMDN(2) 

R*8 

1.0D-5 

Integration lower bounds for gimble angles and 
angular rates 


DMDN(l) (rad) 
DMDN(2) (rad/sec) 

Plotting locations for two axis gimble motion variables 


Input 

Symbol Description 

KPLOTS(3 1 0) A Timuth angle for the gimble azimuth platform 

KPLOTS(31 1) Azimuth angular rate for the gimble azimuth 

platform 


KPLOTS(3 1 2) 


Elevation angle for the gimble elevation platform 


KPLOTS(313) 


Elevation angular rate for the gimble elevation 
platform 
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CONTROL SYSTEM FOR TWO AXIS PLATFORM (DAMPER) MAGNETIC TRACKING 


Fortran 

Symbol Type 
1PLDCSO-20) 1*4 

IPLDCS(I) 


IPLDCS(2) 


IPLDCS(3) 

IPLDCS(4) 

IPLDCS(5-20) 
DCSPRM(IOO) R*8 
DCSPRM(l-9) 
DCSPRM(’O) 
DCSPRM(II) 

DCSPRM( 1 2) 


Preset 

Value Description 

0 Control words for gimble damper platform control 

system simulation 

Control word to invoke the two axis platform control 
system option 

IPLDCS(l) 0 No control system 
IPLDCS(l) 1 Control system 

NOTE: If a control system is invoked the gimble 
simulation must also be invoked 

Control word to specify magnetometer is mounted 
at the tip of an element 

IPLDCS(2) 0 Magnetometer in body 
IPLDCS(2) K Magnetometer on Kth element 

NOTE : K is defined in terms of the internal program 
element ordering system 

Control word for guassian noise transfer function for 
the three axes of the magnetometer 

IPLDCS(3) 0 No transfer function 
IPLDCS(3) 1 Transfer function 

Number of magnetometer samples to be averaged for 
smoothing of magnetometer output 

1PLDCS(4) N average N samples 
N > 1 

Not used 

0.0D0 Platform control system parameters 
Not used 

Sampling time lor PID digital controlled (see) 

Azimuth axis quantatisation level for platform 
position (rad) 

Azimuth PID integrator upper saturation limit 
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CONTROL SYSTEM FOR TWO AXIS PLATFORM (DAMPER) MAGNETIC TRACKING 


Fortran 

Symbol Type 

Preset 

Value 

Description 

DCSPRM(13) R*8 

0.0D0 

Azimuth PL) integrator lower saturation limit 

DCSPRM(14) 


Azimuth PID proportional gain K t > 

DCSPRM(15) 


Azimuth PID integrator gain K[ 

DCSPRM(16) 


Azimuth PID derivative gain Kp 

DCSPRM( 17-20) 


Not used 

DCSPRM(21) 


Elevation axis quantatisation level for platform 
position (rad) 

DCSPRM(22) 


Elevation PID integrator upper saturation limit 

DCSPRM(23) 


Elevation PID integrator lower saturation limit 

DCSPRM(24) 


Elevation PID proportional gain Kp 

DCSPRM(25) 


Elevation PID integrator gain Kj 

DCSPRM(26) 


Elevation PID derivative gain Kq 

DCSPRM(27-29) 


Not used 

DCSPRM(30) 


Integration upper bound for magnetometer first 
order lag transfer function 

DCSPRM(31) 


Integration lower bound for magnetometer first 
order lag tiansfer function 

DCSPRM(32) 


Magnetometer one axis bandwidth 

DCSPRM(33) R*8 

0.0D0 

Magnetometer two axis bandwidth 

DCSPRM(34) 


Magnetometer three axis bandwidth 

DCSPRM(35) 


Magnetometer sampling rate (sec) 

DCSPRM(36) 


Computational delay (sec) 

DCSPRM(37-40) 


Not used 
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CONTROL SYSTEM FOR TWO AXIS PLATFORM (PAMPER) MAGNETIC TRACKING (CONT) 


Fortran 

Symbol 

DCSPRMC41) 

DCSPRM(42) 

DCSPRM(43) 

DCSPRM(44) 

DCSPRM(45) 

DCSPRM(46) 

DCSPRM(47) 

DCSPRM(48-! 

DCSPRM(51) 

DCSPRM(52) 

DCSPRM(53) 

>:SPRM(54) 

DCSPRM(55) 

DCSPRM(56) 

DCSPRM(57) 

DCSPRM(58- 

DCSPRM(80) 


Preset 

Type Value Description 

R*8 0.0D0 Azimuth platform amplifier gain K A 

Azimuth platform motor torque constant K-p 

Azimuth platform back EMF constant Kg 

Azimuth platform motor torque upper limit Lj 
(ft-lbs) 

Azimuth platform motor torque lower limit Li 
(ft-lbs) 

Azimuth platform coulomb friction torque constant 
(ft-:bs) 

Azimuth platform minimum angular rate fcr coulomb 
friction torque (rad/sec) 

r 0) Not used 

Elevation platform amplifier gain K A 

Elevation platform motor torque constant Kx 

Elevation platform back EMF constant Kg 

Elevation platform motor torque upper limit Lj 
(ft-lbs) 

Elevation platform motor torque lower limit L| 
(ft-lbs) 

Elevation platform coulomb friction torque constant 
(ft-lbs) 

R*8 0.0D0 Elevation platform minimum angular rate for 

coulomb friction torque (rad/sec) 

79) Not used 

Amplitude of sinusoidal noise added to one axis 
magnetometer measurement (vauss) 

Amplitude of sinusoidal noise for two axis (gauss) 
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CONTROL SYSTEM FOR TWO AXIS PLATFORM (DAMPER) MAGNETIC TRACKING (CONT) 


Fortran 

Symbol Type 

Preset 

Value 

Description 

DCSPRM(82) 


Amplitude of sinusoidal noise for three axis (gruss) 

rCSPRM(83) 


Phase of sinusoidal noise added to one axis magneto- 
meter measurement (deg) 

DCSPRM(84 ) 


Phase of sinusoidal noise for two axis (deg) 

C°SPRM(85) 


Phase of sinusoidal noise for three axis (deg) 

DCSPRM(86; 


Frequency of sinusoidal noise added to one axis 
magnetometer measurement (cps) 

DCSF’RM(87) 


Frequency of sinusoidal noise for two axis (so... 

DCSPRM(88) 


Frequency of sinusoidal noise for three axis (cps; 

DCSPRM(89) 


Not used 

DCSPRM(90-95) 


Used internally. Not input 

DCSPRM(96-100) 


Not used 

Plotting locations for damper Platform control system variables 

Input 

Symbol 


Description 

KPLQTS(314) 


Output of first order lag transfer function for mag- 
netometer one axis 

KPLOTS(315) 


Output of first order lag 'ransfer function for magnet- 
ometer two axis 

KPLOTS(316) 


Output of first order lag transfer func tion f cr magnet- 
ometer three axis 

KPLOTS(317) 


Azimuth error output 

KPLOTS(318) 


Elevation error output 

KPLOTS(319) 


Azimuth PID digital controlle~ output 

KPLOTS(320) 


Elevation PID digital contro! 1 .* output 

KPLOTS(321) 


Azimuth platform drive motor torque 
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CONTROL SYSTEM FOR TWO AXIS PLATFORM (DAMPER) MAGNETIC TRACKING (CONT) 


Input 

Symbol 

KPLOTS(322) 


Description 


Elevation platform drive motor torque 


(LEFT BLANK) 
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AXIS MOMENTUM WHEEL CONTROL 


Fortran 

Symbol Type 
IAMWHO-10) 1*4 
IAMWH(l) 

IAMWH(2) 

IAMWH(3) 

IAMWH(4) 

AMWHPR (200) R*8 

AMWHPR(l-4) 

AMWHPR(5) 

AMWHPR(b-ll) 

AMWHPR(12) 

AMWHPR(13) 

AMWHPR(14) 

AMWHPR(15) R*8 

AMWHPR(16) 

AMWHPR(17) 


Preset 

Value Description 

0 Control words for axis momentum wheel simulation 

Basic control word for axis momentum wheel 
simulation 

IAMWH(1) 0 No simulation 
IAMWH(l) 1 Simulation 

Control word for roll axis momentum wheel 
1AMWH(2) 0 No roll axis wheel 
1AMWH(2) 1 Roll axis wheel 

Control word for pitch axis momentum wheel 
IAMWH(3) 0 No pitch axis wheel 
IAMWH(3) 1 Pitch axis wheel 

Control word for Yaw axis momentum wheel 
IAMWF'4) 0 No Yaw axis wheel 
IAMWH(4) 1 Yaw axis wheel 

0.0D0 Physical parameters for axis momentum wheel 
simulation 

Not used 

Expenential decay coefficient for integrator satur- 
ation simulation 

Not used 

Upper saturation limit for roll axis integrator 
Lower saturation limit for roll axis integrator 
Proportional gain in roll axis controller 
O.ODO Integrator gain in roll axis controller 
Derivative gain in roll axis controller 
Filter gain for roll angle error 
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AXIS MOMENTUM WHEEL CONTROL (CONT) 


Fortran 

Symbol Type 

Preset 

Value 

Description 

AMWHPR08) R*8 

0.0D0 

Filter bandwidth for roll angle error 

AMWHPRf 1 9-2 1 ) 


Not used 

AMWHPRC2) 


Upper saturation limit for pitch axis integrator 

AMWHPRC3) 


Lower saturation limit for pitch axis integrator 

AMWHPRC4) 


Proportional gain in pitch axis a ntroller 

AMWHPRC5) 


Integrator gain in pitch axis controller 

AMWHPRC6) 


Derivative gain in pitch axis controDer 

AMWHPRC7) 


Filter gain for pitch angle error 

AMWHPRC8) 


Filter bandwidth for pitch angle error 

AMWHPR( 29-31) 


Not used 

AMWHPR(32) 


Upper saturation limit for yaw axis integrator 

AMWHPR(33) 


Lower saturation limit for yaw axis integrator 

AMWHPR(34) 


Proportional gain in yaw axis controller 

AMWHPR(35) 


Integrator gain in yaw axis controller 

AMWHPR136) 


Derivative gain in yaw axis controller 

AMWHPR(37) R*8 

0.0D0 

Filter gain for yaw angle error 

AMWHPR(38) 


Filter bandwidth for yaw angle error 

AMWHPR( 39-40) 


Not used 

AMWHPR(41) 


Roll axis amplifier gain for controller output 

AMWHPR(42) 


Roll axis momentum wheel motor torque constant 

AMWHPR(43) 


Roll axis momentum wheel motor back EMF 
constant 
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AXIS MOMENTUM WHEEL CONTROL (PONT) 


Fortran 

Symbol Type 

Preset 

Value 

Description 

AMWHPR(44) R*8 

0.0D0 

Roll axis momentum wheel torque upper limit 

AMWHPR(45) 


Roll axis momentum wheel motor torque lower 
limit 

AMWHPR(46) 


Roll axis momentum wheel coulomb friction torque 
constnat 

AMWHPR(47) 


RoD axis momentum wheel minimum angular rate 
for coulomb friction torque 

AMWHPR(48) 


Roll axis momentum whce' moment of inertia 

AMWHPR( 49-50) 


Not used 

AMWHPR(5 1 ) 


Pitch axis amplifier gain for controller output 

AMWHPR(52) 


Pitch axis momentum wheel motor torque constant 

AMWHPR(53) 


Pitch axis momentum wheel motor back EMF 
constant 

AMWHPR(54) 


Pitch axis momentum wtieel motor torque upper 
limit 

AMWHPR(55) 


Pitch axis momentum wheel motor torque lower limit 

AMWHPR(56) R*8 

0.0D0 

Pitch axis momentum wheel coulomb friction torque 
constant 

AMWHPR(57) 


Pitch axis momentum wheel minimum angular rate 
for coulomb friction torque 

AMWHPR(58) 


Pitch axis momentum wheel moment of inertia 

AMWHPR(59-60) 


Not used 

AMWHPR(61) 


Yaw axis amplifier gain for controller output 

AMWHPR(62) 


Yaw axis momentum wheel motor torque constant 

AMWHPR(63) 


Yaw axis momentum wheel motor back EMF 
constant 
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AXIS MOMENTUM WHEEL CONTROL (CONT) 


Fortran 

Symbol Type 

Preset 

Value 

Description 

AMWHPR(64) R*8 

O.ODO 

Yaw axis momentum wheel motor torque upper 
limit 

AMWHPR(65) 


Yaw axis momentum wheel motor torque lower 
limit 

AMWHPR(66) 


Yaw axis momentum wheel coulomb friction torque 
constant 

AMWHF 1(67) 


Yaw axis momentum wheel minimum angular rate 
for coulomb friction torque 

AMWHPR(68) 


Yaw axis momentum wheel moment uf inertia 

AMWHPR(69-I00) 


Not used 

AMWHPR(lOl) 


Integration upper bound for filtered angle error 

AMWHPR(102) 


Integration lower bound for filtered angle error 

AMWHPR(103) 


Integration upper bound for momentum wheel speed 

AMWHPR(104) 


Integration lower bound for momentum wheel 
speed 

AMWHPR( 105-1 10) 


Not used 

AMWHPR( 111) 


Initial condition for roll axis momentum wheel speed 

AMWHPR( 112-1 20) 


Not used 

AMWHPR(121) 


Initial condition for pitch axis momentum wheel 
speed 

AMWHPR(122-I30) 


Not used 

AMWHPR(131) 


Initial condition for yaw axis momentum wheel 
speed 

AMWHPR( 132-200) 


Not used 
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AXIS MOMENTUM WHEEL CONTROL (CONT) 


Plotting locations for axis momentum wheel control 


Input 

Symbol Description 


KPLOTS(323) 

Roll axis filler output of error signal 

KPLOTS(324) 

Pitch axis filter output of error signal 

KPLOTS(325) 

Yaw axis filter output of error signal 

KPLOTS(326) 

Roll axis momentum wheel speed 

KPLOTS(327) 

Pitch axis momentum wheel speed 

KPLOTS(328) 

Yaw axis momentum wheel speed 


(LEFT BLANK) 
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ARBITRARILY ORIENTED MOMENTUM WHEEL CONTROL 


Fortran 

Symbol Type 

Preset 

Value 

Description 

IAMPRM(IO) 1*4 

0 

Control words for arbitrarily oriented momentum 
wheel simulation 

IAMPRM(l) 


Control word to invoke simulation 
IAMPRM( 1 ) 0 No simulation 
IAMPRM(l) l Simulation 

AMPARM(200) R*8 

0.0D0 

Physical constants for arbitrarily oriented momentum 
wheel simulation 

AMPARM(l) 


Not used 

AMPARM(2-4) 


Direction cosines of the momentum wheel axis in 
body frame 

AMPARM(5-7) 


Weighting factors for roll(5), Pitch(6) and YAW(7) 
angles input to filter 

AMPARM(8) 


Exponential decay coefficient for integrator satura- 
tion simulation 

AMPARM(9-31) 


Not used 

AMPARM(32) 


Upper saturation limit for integrator 

AMPARM(33) 


Lower saturation limit for integrator 

AMPARM(34) 


Proportional gain in controller 

AMPARM(35) 


Integrator gain in controller 

AMPARM(36) 


Derivative gain in controller 

AMPARM(37) 


Filter gain for angle error 

AMPARM (38) 


Filter bandwidth for angle error 

AMPARM(39-60) 


Not used 

AMPARM(61) 


Amplifier gain for controller output 

AMPARM(62) 


Momentum wheel motor torque constant 
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ARBITRARILY ORIENTED MOMENTUM WHEEL CONTROL (CONT) 


Fortran 

Symbol Type 
AMPARM(63) R*8 
AMPARM (64) 
AMPARM(65) 
AMPARM(66) 
AMPARM(67) 

AMPARM(68) 

AMPARM(69-1 00) 

AMPARM(lOl) 

AMPARM( 1 02) 

AMPARM(103) 

AMPARM( 1 04) 

AMPARM(105-130) 

AMPARM(131) 

AMPARM( 132-200) 

Plotting locations for 

Input 

Symbol 

KPLOTS(329) 

KPLOTS(330) 


Preset 

Value 


0.0D0 


rbitrarily 


Description 

Momentum wheel motor back EMF constant 

Momentum wheel motor torque upper limit 

Momentum wheel motor torque lower limit 

Momentum wheel coulomb friction torque constant 

Momentum wheel minimum angular rate for coulomb 
friction torque 

Momentum wheel moment of inertia 
Not used 

Integration upper bound for filtered angle error 
Integration lower bound for filtered angle error 
Integration upper bound for momentum wheel speed 
Integration lower bound for momentum wheel speed 
Not used 

Initial condition for momentum wheel speed 
Not used 

oriented momentum wheel 


Description 

Filter output of control system angular error 
Momentum wheel speed 
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END OF INPUT 


The end of the input cards is indicated by two cards: the first one with a 1 
punched in column one; the second one with the word END punched in columns 
one through three. 


100 



PART 4 

DEFINITION OF PRINTED COMPUTER OUTPUT DATA 


Fortran 


Symbol 

Description 

Units 

DATE 

Year, month, day 

YY MM DD 

TIME 

Time, in hours, minutes, seconds 
from start of day. 

HHMMSS.XX 

XSAT1 

XSAT2 

XSAT3 

Aries inertial components of satellite 
position vector. 

km 

XSATDT1 

XSATDT2 

XSATDT3 

Aries inertial components of satellite 
velocity vector. 

km/ sec 

DELTAT 

Integration step size at the time of print out 

sec 

SA(i,j) 

Transformation matrix from body frame 
to Aries frame. 

N.D. 

RMAG 

Magnitude of satellite position vector . 

km 

LAT 

Geodetic latitude 

deg 

LONG 

Longitude, measured east from 
Greenwich Meridian 

deg 

ALFAE 

BETAE 

GAMAE 

Euler angles; orientation of gravity 
gradient satellite body axes with re- 
spect to local vertical frame. Output 
only for earth-oriented satellite. 

deg 

PSI1 

PHI1 

THET1 

Euler angles; orientation of spin stabi- 
lized satellite body axes with respect to 
local inertial frame . Output only for 
spin stabilized satellites. 

deg 
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Fortran 

Symbol 

W1BC 

W2BC 

W3BC 


RRATE 

PRATE 

YRATE 

W1B 

W2B 

W3B 

DOUT(k, j) 


DOUTDT(k, j) 


DIN(k, j) 


DINDOT(k, j) 


A(k, j) 


Description 


Body frame components of satellite 
angular velocity measured with re- 
spect to the local vertical frame. Out- 
put only for earth-oriented satellite. 

Roll, pitch, and yaw angle rates. Output 
only for earth-oriented satellite. 


Body frame components of satellite angu- 
lar velocity vector measured with re- 
spect to inertial space. 

Component of the element (libralion 
damper) tip displacement in the j th 
bending modes as measured in the ele- 
ment frame (Z 1 ) along the 2 axis of the 
frame. (INOPT = 2 only) 

Component of the k^ element (libration 
damper) tip velocity in the j 111 bending 
mode as measured in the element frame 
(Z 1 ) along the 2 axis of the frame. 

(INOPT = 2 only) 

Component of the k^ element (libration 
damper) tip displacement in the j 111 
bonding mode as measured in the ele- 
ment frame (Z 1 ) along the 3 axis of the 
frame. (INOPT = 2 only) 

Component of the k 1 * 1 element (libration 
damper) tip velocity in the j 1 * 1 bending 
mode as measured in the element frame 
(Z 1 ) along the 3 axis of the frame. 
(INOPT = 2 only) 

Component of the k* 11 primary element 
tip deflection in the j th bending mode as 
measured along the 2 axis of the element 
frame. 


Units 
deg/ sec 

deg/ sec 
deg/ sec 
feet 

ft/sec 

feet 

feet 

feet 
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Fortran 

Symbol 

Description 

Units 

ADOT(k, j) 

Component of the k 1 * 1 primary element 
tip velocity in the j 1 * 1 bending mode as 
measured in the element frame along 
the 2 axis of the frame . 

ft/sec 

B(k, j) 

Component of the k 1 * 1 primary element 
tip deflection in the j* bending mode as 
measured along the 3 axis of the ele- 
ment frame. 

i^et 

BDOT(k, j) 

Component of the k A primary element 
tip velocity in the j 1 * 1 bending mode as 
measured along the 3 axis of the ele- 
ment frame. 

ft/sec 

UD(i, k) 

Total tip displacement, for the bending 
modes simulated, of the k** 1 libration 
damper element , measured along the 
i 1 * 1 axis of the element frame. 

(INOPT = 2 only) 

feet 

UD(i)DT(k) 

Total tip velocity, for the bending modes 
simulated, of the k 1 * 1 libration damper 
element, measured along the i 1 * 1 axis of 
the element frame. (INOPT = 2 only) 

ft/sec 

U(i, k) 

Total tip displacement, for the bending 
modes simulated, of the k^ primary 
element, measured along the 1 th axis of 
the element frame. 

feet 

U(i)DOT(k) 

Total tip velocity, for the bending modes 
simulated, of the k* primary element, 
measured along the I th axis of the ele- 
ment frame . 

ft/sec 

ZLD(k) 

Length of the k* element, libration 
damper. (INOPT = 2 only) 

feet 

ZLK(k) 

Length of the k 111 primary element. 

feet 


103 



Fortran 

Symbol 

PHILD 


DPHILD 


OMEGL 


VSUBL 


MSUBM1 

MSUBM2 

MSUBM3 

CMX 


SIMPX 


SUNVEC1 

SUNVEC2 

SUNVEC3 

XMB1 

XMB2 

XMB3 

SMAGI1 

SMAGI2 

SMAGI3 


Description 

The angular deflection of the libration 
damper boom relative to its equilibrium 
position. (INOPT = 2 only) 

Angular velocity of li nation damper 
boom relative to the body. 

(INOPT = 2 only) 

Spatial average angular velocity of the 
nutation damper fluid relative to the 
body. (INOPT = 1 only) 

Average linear velocity of the nutation 
damper fluid relative to the damper tube 
wall. (INOPT = 1 only; 

Body frame components of torque exerted 
by the nutation damper fluid upon the 
satellite. (INOPT = 1 only) 

Control moment along Yj body axis. 
(INOPT = 1 only) 

A function related to total impulse for 
control moments applied along Yi body 
integral of the applied moment. 

(INOI T = 1 only) 

Components of sun to satellite unit 
vector , body frame . 


Body frame components of magnetic 
torque acting on satellite. 


Components of Earth's magnetic field 
strength, Aries inertial frame. 


Units 

deg 

deg/ sec 
deg/ sec 

ft/sec 

ft- lb 

ft- lb 
ft-lb-sec 

N.D. 

ft-lb 

Gauss 


104 



Fortran 

Symbol Description Units 

SMAGB1 Components of the Earth’s magnetic field Gauss 

SMAGB2 strength, body frame. 

SMAGB3 

SOLILL Solar illumination N. D. 

0 = Occulted 

1 = Full sunlight 

RWHEEL1 Body frame components of reaction torque ft-lb 

RWHEEL2 exerted by the momentum wheels upon 

RWHEEL3 the satellite. 

EPSERR Angular error between desired and actual deg 

direction of spin axis. (INOPT = 1 only) 

SUNANG Angle between sunline and orbit normal. deg 

SB(BODY) Sun to satellite unit vector in body frame. N. D. 

SD (INERTIAL) Sun to satellite unit vector in inertial N. D. 

coordinates - print out when boom de- 
ployment from sun pulse is used. 

FTAKIN(10, 3) Thermal force in flexible element 1-2 lb 

plane direction before satellite shadowing 
and thermal lag. 

FTAKOT(10,3) Thermal force in flexible element 1-2 lb 

plane direction after satellite shadowing 
and thermal lag. 

FTBKIN(10,3) Thermal force in B direction before lb 

satellite shadowing or thermal lag. 

FTBKOT(10,3) Thermal force in flexible element 1-3 plane lb 

after satellite shadowing or thermal lag. 

TENSN(IO) Element root tension. lb 

YCEMS(i) Displacement components of center of feet 

i=l, 2, 3 mass with respect to original body 

axes due to flexible appendage motion 
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Fortran 

Symbol 

Math 

Symbol Description 

Unit 

BNMTA(IO) 

Element root bending moment in 
flexible element 1-2 plane. 

ft-ib 

BNMTB(IO) 

Element root bending moment in 
flexible element 1-3 plane. 

ft-lb 

MOMENT1 

MOMENT2 

MOMENT3 

Right hand side of Euler equations of 
motion applied to the spacecraft hub with 
gyroscopic terms also included. Used in 
check out- 

ft-lb 

HAMILT 

Relative Hamiltonian of entire system. 

ft-lb 

HVECTR1 

HVECTR2 

HVECTR3 

{hf } Angular momentum vector of the 

spacecraft in inertial frame 

Ib-ft-sec 

HBODY1 

HBODY2 

HBODY3 

{ hjj } Angular momentum vector 

expressed in spacecraft 
body reference frame 

ft-lb-sec 

HMAG 

Magnitude of the angular momentum 
vector 

Ib-ft-sec 

NUTANG 

9 Angle between the angular momen- 

tum vector and the spin axis (3 axis ) 

deg 

EPSERH 

N.A. Angular error between the initial 

angular momentum vector direction 
and actual direction of spin axis 
(1NOPT = 1 only) 

deg 

CW(k, j) 

1L 

Twist angle of the^ flexible 
element with the j twist mode 

deg 

CWD(k, j) 

Tjjgne derivative of twist angle x>f 
k flexible element w*** the j 
twist mode 

deg/sec 

STAG PR 

Stagnation pressure (IDRAG > 1) 

lbs/ft 2 
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Classical Orbital Parameters After Thrusting 

The change in the orbit velocity vector at the time of thrusting will be reflected 
in the change of the values of classical orbital parameters. The standard 
printout will be interrupted and a set of new orbit information will be printed 
immediately after the impulsive thrusting. If the sun reference option is used 
(ISPLSE=1), the number of pulses together with the mean time (hour, minute, 
second) for the orbit update will always be printed. In this case, the orbit up- 
date message printout will be controlled by input word ISPNP. For instance, 
if the IPULSE is specified as 21 and ISPNP=5, the orbit update message will 
be printed at 5, 10, 15, 20, and 21st pulse. 

Fortran Math 

Symbol Symbol Description Units 

ACCOB(I) {r} Acceleration vector of center ft/ sec 2 

of mass expressed in Aries 
inertial frame (I = 1, 2, 3) 

ACCRED(J) a hl Jth accelerometer reading ft/ sec 2 

(J = 1, NUMHUB) 

HUBACC(I.J) Acceleration vector ex- ft/sec 2 

pressed in Aries inertial 
frame of a point in the hub 
(I = 1,2,3; J = 1, NUMHUB) 


TIPACC(I.K) 

{a-j} Acceleration vector of kth 

element tip expressed xn 
Aries inertial frame 
(1 = 1,2,3; K = 1, NELMTS) 

ft/ sec 2 

BIXX /• 

•'s 

(y 1 - + v 2 1 

\ 2 Lp '3L Pi 

Jd in 

Instantaneous moment of 
inertia with respect to 
No. 1 body reference axis 

slug- ft 2 

BIYY r 

*'s 

/y 2 + y 2 1 

\ 1 Lp 3 L pj 

jdm 

Instantaneous moment of 
inertia with respect to 
No. 2 body reference axis 

slug-ft 2 


BIZZ 



Instantaneous moment of 
inertia with respect to 
No. 3 body reference axis 


slug-ft 2 
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Fast Fourier Transform Analysis Output 


The related printer plots for the Fast Fourier Transform are the power spec- 
trum density function and the cross-correlation functions together with the re- 
sults of harmonic analysis (frequencies, amplitude and phase angles). These 
are printed immediately following the regular printer plot of the FSD program. 
The data set to be analyzed and the autocorrelation function plot are also avail 
able. When the Fast Fourier Transform analysis option is used, increasing the 
I/O time estimate by a factor of 50 percent is recommended. 

DE-B Control System Printed Output 


Fortran 

Symbol 

Description 

Units 

PTCH Out 

Pitch channel sensor output 

volts 

ROLL Out 

Roll channel sensor output 

volts 

COMP Out 

Output of compensation amplifier 

volts 

TACH Out 

Output of tachometer 

volts 

TMOTOR 

Torque output of momentum wheel motor 

ft.-lbs 

WHL SPD 

Momentum wheel speed 

rad /sec 

NUTD Out 

Nutation damper phase shift circuit output 

olts 

Secondary Body 

Printed Output 


GAMSB 

Rotation of secondary body about 3 axis 

deg 

ALPSB 

Rotation of secondary body about the carried 
1 axis 

deg 

BETSB 

Rotation of secondary body about the carried 
2 axis 

deg 

OM1SB 

Component of secondary body relative angular 
velocity on 1 axis of secondary body 

deg/sec 

OM2SB 

Component of secondary body relative angular 
velocity on 2 axis of secondary body 

deg/sec 

OM3SB 

Component of secondary body relative angular 
velocity on 3 axis of secondary body 
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Fortran 

Symbol 


Description 


Units 


SEN OUTP Sensor outpu' (zero crossing of state variable) depends on 

state 

FIL OUTP Filter output (zero crossing of state variable) variable 


TWO AXIS PLATFORM (and DAMPER) CONTROL SYSTEMS PRINTED OUTPUT 


Fortran 
S> nbol 

Description 

Units 

Stf iGFl 
SMAGF2 
SMAGF3 

Filtered body frame components of earth’s 
magnetic field 

Gauss 

GMBL AZ 

Gimble platform azimuth angle 

d( ? 

GMBL AZD 

Gimble platform azimuth angular velocity 

deg/sec 

GMBL EL 

Gimble platform elevation angle 

deg 

GMBL ELD 

Gimble platform elevation angular velocity 

deg/sec 

AZIM ERR 

Gimble platform azimuth error angle 

deg 

ELEV ERR 

Gimble platform elevation error angle 

deg 

AZIM PID 
ELEV PID 

Gimble platform azimuth PID output 
Gimble platform elevation PID output 

Depends on PID 
gain constants 

AZIM MOT 

Gimble platform azimuth motor torque 

ft-lbs 

ELEV MOT 

Gimble platform elevation motor torque 

ft-lbs 

AXIS MOr 

N TUM WHEEL PRINTED OUTPUT 


For an 
f mbol 

Description 

Units 

FTR LOLL 

Roll filtered error output 

deg 

FTR PITCH 

Pitch filtered error output 

deg 

FTR \ AW 

Yaw filtered error output 

deg 
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\RBITRARILY ORIENTED MOMENTUM WHEEL PRINTED OUTPUT 


Fortran 


Symbol 

Description 

Units 

FTRANG 

Filtered error output 

deg 

ARB MWS 

Momentum wheel speed 

deg/sec 

MOTTK 

Motor torque 

ft-lbs 

CONT INT 

Output of PID control system integrator 

Depends on PID 
gain constants 

CONT DER 

Output of PID control system differentiator 


ELEMENT THERMAL EXPANSION PRINTED OUTPUT 


Fortran 

Symbol 

Description 

Units 

TEMP(1-10) 

Average temperature of element 

°R 

ZLKPO-10) 

Axial velocity of element due to temperature effects 

ft/sec 

ZLKDPC1-10) 

Axial acceleration of element due to temperature 
effects 

ft/sec 2 
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Adams-Moulton Numerical Integration Control Message Output 


The Adams-Moulton numerical integration control message is presented in a 
table format immediately following the nominal FSD printout and before (be 
printer plot in each stacked case sequence. Ibis table is the summary of the 
performance of numerical integration within a stacked case run. The definitions 
of the variables in the table are as follows: 


Variable 

Description 

KCUT 

Number of times the integration step size is 
reduced for a particular equation 

UPPER 

Upper integration bound for a particular equation 

LOWER 

Lower integration bound for a particular equation 

DEPEND 

Dependent array value of a state variable at the 
starting time of a stacked case sequence 

DERIVE 

Derivative array value of a state variable at the 
starting time of a stacked case sequence 


111 




GENERAL REFERENCES 


1. Avco Corporation, AVSD-0191-71-CR, User's Manual for IMP Dynamics 
Computer Program, Volume I, E. A. Lawlor, L. Beltracchi, L. Turner, 
and M. Weinberger, March 1971. 

2. Avco Corporation, AVSD-0191-71-CR, User's Manual for IMP Dynamics 
Computer Program, Volume n. Integral Evaluation Computer Program, 

E. A. Lawlor, March 1971. 

3. Avco Corporation, Contract No. NAS 5-24008 Mod 15, Modification of the 
FSD Program for Linear Thrusting and Angular Momentum Computation 
User's Manual , K. Yong and E. A. Lawlor, November 1974. 

4. Avco Corporation, Spacecraft Configuration Plot Program, A. Anderson, 
February 1973. 

5. X-732-73-151, A User's Guide to the Flexible Spacecraft Dynamics 
Program- 1 , August 1973. 

6. E. G. Stassinopolus, G. D. Mead, ALLMAG, GDALMG, L1NTRA, Computer 

Programs For Geomagnetic Fields and Field line Calculations, NSSDC -72-12 
FEB 1972. ” 

7. X-7 12-76-4, A User's Guide to the Flexible Spacecraft Dynamics 
Program-n, March 1976 

8. Avco Corporation, Modification of the IMP Dynamics Computer Program for Limited 
Dual Spin Capability, E.A. Lawlor, February 1972. 


112 




APPENDIX A 

COORDINATE SYSTEMS AND 
OTHER RELATED MATERIAL 



DIRECTION: 
NORTH POLE 



ECLIPTIC PLANE T 

(APPARENT SUN ORBIT) " 


DIRECTION: 
VERNAL EQUINOX 


ARIES FRAME: 

RIGHT HANDED ORTHOGONAL REFERENCE FRAME 
ORIGIN: CENTER OF EARTH 

I„ POINTS TOWARD VERNAL EQUINOX 

I 2 , EQUATORIAL PLANE, PERPENDICULAR TO T, L 

T 3 , POINTS TOWARD EARTH'S GEOMETRIC NORTH POLE 


THE LOCAL INERTlAL (L-FRAME) IS OBTAINED BY A PARALLEL TRANSLATION 
OF THE ARIES FRAME FROM THE EARTH CENTER TO THE SATELLITE CENTER 
OF MASS. 

Figure 1. Aries Coordinate System 
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L: LOCAL INERTIAL FRAME 
Y: BODY FRAME 

SEQUENCE OF ROTATIONS 4 , ^ e -* <p 

'.LONG AXIS (3) (1) (3) 

FOR SPINNING BODY: Y 3 TAKEN AS SPIN AXIS 

0 * ANGLE BETWEEN SPIN AXIS AND THIRD INERTIAL DIRECTION 

Figure 2, Euler Angles for a Spinning Body 
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PERPENDICULAR TO 
INSTANTANEOUS 
ORBITAL PLANE 



lx POINTS FROM EARTH CENTER TO CENTER OF MASS OF SATELLITE, 
AND IS PARALLEL TO THE RADIUS VECTOR I 

I 2C PERPENDICULAR TO INSTANTANEOUS QRBITAL PLANE 

ITS DIRECTION IS DEFINED BY x V. V IS THE INERTIAL VELOCITY 
OF THE CENTER OF MASS OF THE SATELLITE. 

T, c COMPLETES RIGHTHANDED ORTHOGONAL FRAME 
FOR A CIRCULAR ORBIT T,c LIES ALONG V 
FOR AN ELLIPTIC ORBIT, I 1C LIES ALONG V ONLY AT APOGEE 
AND PERIGEE. T, c ALWAYS LIES IN INSTANTANEOUS ORBITAL 
PLANE AND MAKES AN ACUTE ANGLE WITH \7 

THE LOCAL "VERTICAL" FRAME (C) HAS UNIT VECTORS T 1C , I 2C , Tac AND 
IS CENTERED AT CENTER OF MASS OF SATELLITE 

THE "VERTICAL" USED HERE IS DIRECTION OF GRAVITY FORCE FOR A 
SPHERICAL EARTH. 


Figure 3. Lc al Vertical Coordinate System 
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*£ ’ n ROTATION 


DIRECTION COSINES OF D WITH RESPECT TO LOCAL 
INERTIAL FRAME: 

SIN tj COS f 
SIN r? SIN f 
COSr? 


Figure 4. Error Angle for Spin Axis 





Y 3 



yi, Y2, y 3 COMPONENTS OF AN ARBITRARY VECTOR, ON THE 
>, 2, 3 AXES OF THE BODY Y FRAME 

z u z 2 , z 3 CC.V.T ONENTS OF AN ARBITRARY VECTOR, ON THE 
1, 2, 3 AXES OF THE LIBRATION DAMPER Z FRAME 


Figure 6. Origin of Libration Damper Z Frame 
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Y 3 



\ 


y 1r y 2 , y 3 COMPONENTS OF AN ARBITRARY VECTOR, ON THE 

1, 2, 3 AXES OF THE BODY Y FRAME 

z .« , z„ , z,, COMPONENTS OF AN ARBITRARY VECTOR, ON THE 
1, 2, 3 AXES OF THE LIBRATION DAMPER Z FRAME. 
Z„ POINTS IN AXIAL DIRECTION OF DAMPER 
BOOM, OUTWARD FROM SATELLITE CORE 

■*1.0 Ad ' 0ld EULER ANGLES 


Figure 7. Orientation of the Libration Damper Z Frame 
with Reppect to Body Y Frame 
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y,, y 2 , y 3 COMPONENTS OF AN ARBITRARY VECTOR, ON THE 

1, 2, 3 AXES OF THE Y FRAME 

z’ 1K , z 3 k COMPONENTS OF AN ARBITRARY VECTOR, ON THE 

I, 2, 3 AXES OF THE Z\ ik* ELEMENT) FRAME. THE 
1 AXIS OF THE ELEMENT FRAME POINTS ALONG UN- 
DEFORMED ELEMENT, OUTWARD FROM SATELLITE 
CORE 


yiOK , y 2 ok , yaoK Y FRAME COMPONENTS OF Z ORIGIN 


Figure 8. Origin of Boom Element Frame 
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Z',k 


y,, y 2 , y 3 COMPONENTS OF AN ARBITRARY VECTOR, ON THE 

I, 2, 3 AXES OF THE REFERENCE Y FRAME 

z* tk , z* 2k , z' 3k COMPONENTS OF AN ARBITRARY VECTOR. ON THE 
1, 2, 3 AXES OF THE Z' k (K th ELEMENT) FRAME. 

THE 1 AXIS OF THE ELEMENT FRAME POINTS ALONG 
UNDEFORMED ELEMENT, OUTWARD FROM SATELLITE 
CORE 

o k , fi k , 7 k EULER ANGLES 

Figure 9. Orientation of Element Frame (Z*) with Respect 
to Reference Y Frame 
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orig<mal - lj 

OF POOR QUALITY 



f,„, z ?k , z Jk COMPONENTS OF AN ARBITRARY VECTOR, ON THE 
1, 2, 3 AXES OF THE Z\ (W"' ELEMENT) FRAME 

u j ELEMENT TIP DEFLECTION, Z Jk COMPONENT, 

u 3 ELEMENT TIP DEFLECTION, Z' 3k COMPONENT 

Figure 10. k tfl Element Coordinate System 


ELEMENT 



Figure 11 


ELEMENT 
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ORIGINAL : s . 
OF POOR QUALITY 


z 



X, Y, Z SATELLITE BODY AXES 

DIRECTION COSINES OF ^ WITH RESPECT TO BODY 
T? x = SIN V D COS?D 
r? v = SIN i)o SINfo 
= COStj 0 

Figure 12. Axis of Ring Damper rj Unify Vector Perpendicular 
to Plane of Ring Damper 


A-ll 


CfJ.TT!."’; 

OF FOG.': c- 


Y, 



s,v = o 
S,y < o 


Figure 13. Sun Sensor Reference Point in Spin Period 



T ? = T, + B r 3 
T 3 = T, + C t 3 

t , = SPIN PERIOD = J-L 

|w M | 

T, = REFERENCE TIME IN SPIN PERIOD 

Figure 14. Time History of Control Moment about Y, 

Body Axis 
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0R:3*:‘ . 

OF POOR £<-•/ 


V 3 



B 1Y/ B 2y , b 3Y components of b y on body axes 

m )Y , m 2Y , m 3Y COMPONENTS OF m Y ON BODY AXES 



DIPOLE MOMENT 
OF SATELLITE 


Figure 15. Dipole Moment of Satellite and Earth's Magnetic Field Flux 

in Body Frame 
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^1L» ^ 2L» ^3L 

UNIT VECTORS, ARIES EQUATORIAL INERTIAL FRAME 

BV: 

VERTICAL COMPONENT OF FIELD, POSITIVE AWAY 
FROM CENTER OF EARTH 

BE: 

PARALLEL TO EQUATORIAL PLANE, TANGENT TO 
CONSTANT LATITUDE CIRCLE, POSITIVE IN EASTERLY 
DIRECTION 

BN: 

COMPLETES RIGHTHANDED ORTHOGONAL TRIAD, 
POSITIVE TOWARDS GEOGRAPHIC NORTH 


Figure 16. Components of Magnetic Field 
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GEOGRAPHIC 

NORTH 



DIRECTION OF 
INTERSECTION 
OF GREENWICH 
MERIDIAN WITH 
EQUATORIAL 
PLANE 


<t> = GEOGRAPHIC LONGITUDE 
6 = GEOGRAPHIC COLATITUDE = 90 6 -LATITUDE 
r = RADIUS VECTOR 


Figure 17. Coordinates for Spherical Harmonics 


. EASTERLY 
DIRECTION 
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Thrust Loading 


The time history of the thrust loading is as shown in Figure 18. 



TIME (Seconds) 


Figure 18. Thrust Time History 


Hie analytic expressions for the thrust loading time history are as follows. 


’ t <ti 

f.ati t, 
t, * t * <« 

1 >t* 


Fity-o 

.eit-t.). 

F.cVbACI-S ] 

F. tr> -- F, (*. ) + 

rt -t 7 

Frt)*o 


a* 


The total impulse under the thrust loading time curve is obtained by integration. 
The impulse is used to apply the Increment in velocity to the trajectory of the 
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spacecraft. The mean time, at which the velocity increment is applied, is ob- 
tained as follows. 


i 



( 2 ) 


The appropriate expressions for the three pieces of the thrust loading time curve 
are given below. 

For 


i, 


I, 

M, 


si S it 


3 / f. d ~e (1-2 >7 

Jt, 


F,rt.) 

B* 


(3) 


For 

t a * t * f i 





r 




r*> 


(4) 


r F, rf*)<*4* C*l, 
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For 


{ * !<, 



where At 21 = t 2 -t, , At 3| = t 3 -t, , At 41 = t 4 -t, etc. 

Thus, the mean time for which the impulsive thrusting 1° applied is given as 


J — Sf * — 

The magnitude of change of velocity at t is given as 

*v- — - — *— 

Ms 

where M s = Mass of the entire spacecraft. 


(6) 


( 7 ) 


OK&tft/lL r» • 

OF POOR QUALITY 
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HORIZONTAL PLANE 



1 CYCLE 



A-20 



A-21 


RELATIVE WHEEL SPEED 



Figure 22. Dynamic Explorer-B Pitch and Active Nutation Damper Block Diagram 













MATHEMATICAL FORMULATION FOR DUAL SPIN CONTROL SYSTEM SIMULATION 

A block diagram of the system to be simulated is given in Figure 22. This block diagram is 
represented as a set of first order ordinary differential equations, which are integrated in paral- 
lel with the equations of motion for the rest of the spacecraft, using the same time step and 
integration algorithm. (The subroutine ADMIMP.) 

For the most part, the control components are linear dynamic systems. For such components, 
the stated transfer functions have been converted to state variable equations using standard 
techniques. This transformation, however, is not always unique. Hence, it is necessary to 
state the exact form utilized in each case. 


In the equations following, the subscripts 1 . 2 etc. .ire used primarily for convenience. How- 
ever. the ordering of variables is the same as in program code. Hence a fourth order model 
with state variables x, . . . x 4 may appear in program code as x 6 . . . x 9 . The actual sub- 
scripts used in program code are given in the section on program inputs. Also, in this section, 
the symbols u and y denote (respectively) the input to and output from the given transfer 
function. In the system simulation, the blocks are coupled together. 

Pitch or roll sensor (2nd order model) 

The transfer function 

(Vs-Q 1 

is represented as 



y=x, 


Pitch or roll sensor (4th order model) 

The transfer function 

. - 

(ii-zSTs-f T 2 s z r 

' S 5 5 
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is represented as 



y=x, 


Nutation Damper Phase Shift Circuit 
The transfer function 



is represented as 



(the minus sign includes the sign inversion corresponding to path A in Figure 22) 


The inversion of sign at various roll angles and the switch between primary and secondary 
damper may be simulated by changing the input parameters. 

Tachometer 


The transfer function 


is represented as 


Kf 

h- si; 




y=x, 
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For the reference control system, the tachometer is significantly faster than the other dynamic 
elements. This causes the tachometer equation to dominate the time step control in numerical 
integration, while having little effect on system performance. Replacing the above transfer 
function with the static operator y = K F U permits a significant reduction in program execu- 
tion time. This alternative model is optionally available in the modified program. 

Pitch Compensation Amplifier 

In unsaturated operation, the transfer function 

K<i(l + sT,) 

O+sTi) 

is represented as 

ii x > - 

y= x,+ /<<,=£ a 

‘z. 

Saturation occurs if 

/ V/ >• 

If this occurs. Y is replaced by Vjj m sign (Y). 
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PROPORTIONAL 
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I FIGURE 24. Generic Azimuth and Elevation PID Digital Controller for Magnetic Field 
' Tracing Control System 






APPENDIX B 


EXAMPLES OF RAE, IMP-I AND DE-B SPACECRAFT INPUT 



ORIGINAL T 
OF POOR QUALi' 


RAE SPACECRAFT INPUT 


INPUT CARDS READ 

DATA*H******«*** DYNAMICS STUDY FOR RAfc-C *A IK K-SOOM DEPLOYMENT •**••••••••• 

DATA*H************»«A******** OI POLES ARE IK ••••»••*•*••*•**•••••••••••••••• 

DATA* AEARTH 1 736.5700 «W0 BS.CCO FLAT *029.3600 2 MU *.9027760*03 

DATA* J2 0.000 J3 0.000 J* O.CCO J22 C.COO 2 J20 0.000 
DATA* TY6R 6.60*0*06 ECLPTC |. « 250*00 
DATA* (TAPE 0 I OR 6 0 

DATA«h******* INCLINATICN a 6* . 0 OEG. ECCEkTRICITY u O.OOS 

DATA* A SAT 2636.570 0.0 0.0 

DATA* XSATOT 0.0 -0*576120 >1.161222 

OAT A* CO 0*0 AREA 10 T036 MUSS I*. 1925 ICEKS 0 

DATA* YEAR 1973.0 Z MONTH 6.0 CAY 17.C 

DATA* IOATE 730*17 

DATA* TINE 0.0 

DATA* TSTQP 9*6.0 

DATA* FREQ 30.0 

DATA* I PLOT 2 

DATA* INOPT 2 

DATA* I0AMP 0 

DATA* IHAMLT I 

DATA* I GRAY t 

DATA* DELTAT O.COI 

DATA* 60T MI 17.155 0.285 0.0 0.265 16.176 0.0 0.0 0.0 21.953 

DATA* 2 MS I*. 1925 

DATA* SCO 10.0036 

DAT A* NELMTS 6 

DATA* NDAMPR 2 

OATA* ZLO l.C 1.0 1.0 1.0 I.C 1.0 l.C 1.0 

DATA* ZLl 0.0 0.0 0.50 C.50 C.53 C. *0 C .0 0.0 

DATA* ZL * O.C 0.0 0.0 0.0 0.0 C .0 0.0 0.0 

OATA* 282(1*02 1.36* -1*36* I.C9T -I .097 1.097 -1.097 0.0 0.0 

DATA* 262(2.1)2 0.371 -C.371 O.C 0.0 C.O C.O 0.0 0.0 

DATA* ZBZI 3.112 -0.227 -0.227 C.*93 0.A93 -0.363 -0.383 0.0 0.0 

DATA* ALFAEK 0.0 0.0 0.0 0.0 0.0 0.0 0.0 C.O 

DATA* BETAEK C.O 0.0 -*0.C -120. C 60.0 120.0 0.0 0.0 

DATA* GAMAEK 15.0 -16?. 0 0.0 0.0 C.O 0.0 C.O 160.0 

DATA* YIZMI 1*112 0.0 

DATA* V IZM( 2*1)2 C*J 

DATA* YlZMC **4)2 -1*6*6 

DATA* MOOESM-21 0 

DATA* MOOES! 3— 6 > 1 

DATA* MOOES! 7- 6 ) 0 

OATA* CmOOLS! 1 -6 ) 1.1E7 

DATA* RTuOC!l-2l 0.2625 

DATA* RTUBC! 3-6) 0.290 

DATA* RTU6E! 7-6) 0*290 

DATA* HTU6C! 1-6 } 0.0022 

DATA* TmERMCCI-6) 6.65E-6 

DATA* TfPMSIl-2) 0.0 

DATA* TIPMSI3-6) 0.00513 

OATA* TIPNS!T-8) 0.0 

DATA* S*0 ( 1-6 ) *.e30-2 

DATA* *(3.1)1 0.&O0 C.OOC O.OCO O.OOC 

DATA* AOOT! 3.1)1 C.000 O.ODO <.000 O.CCO 

DATA* 813*1)1 O.ODO 0.000 C.CCO 0.000 

OATA* 800T (3*1)1 C.000 0.000 C-000 O.CCC 

DATA* D IN! 7.111 0.000 0.000 

OATA* OINDOT (7.1 II O.OOC C.OOC 

DATA* OOUTI7.I11 0.000 0.000 

DATA* DOUTDT (7.1)1 0.000 C.000 

DATA* STMKI 1—6 I 0.0 

DATA* RHOAO-2) A.5652E-* 

OATA* RMOKO-6) *• 3* 76E- * 

DATA* RH0K(7-6) *. 5963E- * 

DATA* POO 0.0 
DATA* COAMP! 1*1-6' 0.0 
OATA* OT0O 0.0 
OATA* AEPO 0*0 
DATA* OETLO 0.0 
0*’A* CARLO 65*0 
DATA* PHILO 0*0 
OATA* DPHlLO 0.0 
OATA* PHIS 35.0 
DATA* ZK10 1 .036E-02 
DATA* ZK20 2.0 
OATA* OECAV 20.0 
OAT AA ZMOO 0.96E-03 
DATA* ZMOBO 0.0 

OATA* OPRM l(l.l)2 0.1 0.0 0.0 

f)A*A* OPRM 1(2.1)? 0.0 0.1 0.0 

DATA* 0PRMKJ.1I2 0.0 0.0 0.1 

DATA* ALFAF 0.0 
DATA* SETAE 30.0 
DATA* GAMAE 196. c 

DATA* OMBCII) 0.0 OMBC ! 2 ) -0.02666 
DATA* OM0CC3) 1.3 

OATA* MLAST 0 MOPLV l OOPLV 0 

oata*i 


• OATA 
•DATA 
•DATA 
•DATA 
•OATA 
•OATA 
•OATA 
•OATA 
•OA TA 
•OATA 
•OATA 
•OATA 
•OATA 
•OATA 
•OATA 
•OATA 
•OATA 
•OATA 
•OATA 
•OATA 

• OATA 
•OATA 
•OATA 

• OATA 

• DATA 
•DATA 
•OATA 
•DATA 
•OATA 
•OATA 
•OATA 

• OATA 

• OATA 
•DATA 

• OATA 
•CATA 

• OATA 
•OATA 
•DATA 
•OATA 

• OATA 

• OATA 

• OATA 
•OATA 

• OATA 

• OATA 
•OATA 
•OATA 
•OATA 
•OATA 
•OATA 

• OATA 
•OATA 

• OATA 
•OATA 
•OATA 
•OATA 
•OATA 
•OATA 
•OATA 
•OATA 
•DATA 
•OATA 
•OATA 

• OATA 
•OATA 

• OATA 
•OATA 
•OATA 

• OATA 
•OATA 

• OATA 
•OATA 

• OATA 

• OATA 
•OATA 
•DATA 
•OATA 

• DATA 

• OATA 

• DATA 
•OATA 
•OATA 

• DATA 

• OATA 
•CATA 

• OATA 
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OF POOR Qo AL.'.V 


IMP-I SPACECRAFT INPUT 


INFLT C A R C S READ 

DATA* X SAT 0.0*444 2SI.30 12 373*0 

CAT AT ASA1CT C.3Q017SE2 T-5t*£7l4 . 0*13364 

OAT X * CC 2*0 AREA 26*6 VMASS 1S.S<50 IOEK< C 

^&4XA*_JtE6«_J S7L.0 i"Cr.TH C<»Y 21*J?391 

DATA# IOPP 0 
CATA* 1 T APE 0 
DATA* TSTCPO 396CC.0 
OAT A * IC AT EL 7SC221 
CATA« TIME 32626.0 

—CATA* — I5.Y QP 2CC. 

DATA* FREQ 1 .0 
CAT A • IPLCT 2 
OATA# INOP T 1 
CAT A t_ lhAKLT 0 
CAT A 4 IGCAV 0 

. CATA* MJCC£S(lr-ll_i — 

DATA* MCPLV | 

CATA# CELT AT 0.01 
data* factor C.5 

DATA# SCYMI 90.0 0*0 0*C 0. C 90. £ C.C 0.0 0.0 117.0 
DATA# SCC 26.6999 

_£AT> f 

CATA# NELMTS e 

CATA* 2L0 200. 200. 20 C* 200* 1C. 1C. 

DATA* 2 l 1 0* 0. 0. C. C.l 0*1 

CATA* 2eZtl.l)2 1 *8094 7 -1.24210 -1.60947 1*2*210 C.O 0.0 
CAT## 262(2*11* I.2421C i.cK94? -1 • £♦ 31 0 - 1 *0C 94 7 C.C C.C 


--11*22006—4^6 rA*£ — 


DA T A * ALF*EK(l-6> 0*0 

CATA* eETAEK C. _C. 0. 0. -9G. 9C. 

OATA* GAMAEK 24.499 12*. *89 214.469 204.469 0*C 0*0 
CATA# EMCOCSU-AJ 1. S^7 
CATA# EMCCLSC5-6I 2.CE7 

-CAT A# fiTUfl£(l-AJ_A.2S- 

DATA* RTU^ECS-C) C * 56 
OATA* MTUee<l-*l 2.0E-J 
DATA* FTU6E<*-0 * . OE - 3 

OATA# THEfiHCCl-tl 6. £32-6 
CATA* T JP*S< l-C I 0.0 

-OAJ-A * -SACUrfil Q-^QA16- . 

CATA* STMM1-4I 0. 

DATA* RKCK11-4J 4.34££2E-4 
OATA* FHCKC5-C) 2.46-J 

data* pcd i.ase-7 

CATA* CTCC 0.0 

-CAXAf „.CCAHgLLtXl 2 Q * 0.3 l„fl-*Q£-l JaUU C 

CATA* CCA*P(<*1)2 0.CC1 0.0C1 C.CC1 C.OOl C.7 0.7 
DATA# CCAMP( 3.1J2 -C.cOOl C.OCI Q,CCl C*001 C.7 C»7 
DATA* AERO 0-0 

DATA! PSil-D.O IfETl -SflUO F*lXl 9C.0 
OATA* C ME C 0.0 O.C 9 t . 


-DATA-# -EJIA-SD-.G- 
OATA* ISPIN3 C 
CATA* XCAUF 0 
OATA* IV1SCS 0 
CATA* 1 AT TOE 0 _ 
CATA* I »m£6l 0 


ClIA t IVGylS 0 


OATA* CPMAG 0.0 0.0 50. C 
CAT/* MACFLD 4 
OATA* MLAST 0 
CATA* 1 


♦ OATA 

♦ OATA 

♦ OATA 
AG-ATA 

♦ OATA 

♦ OATA 

♦ OATA 
♦OATA 

♦ OATA 

•QATA. 

♦ OATA 

♦ OATA 

♦ OATA 

♦ OATA 

♦ OATA 
♦DATA 

♦ DATA 

♦ OATA 

♦ OATA 

♦ DATA 

♦ OATA 
♦OATA 

♦DATA 

♦ DATA 

♦ OATA 
♦OATA 

♦ OATA 
4UATA. 

♦ DATA 
_ *CAT A 

♦ CATA 

♦ OATA 

♦ CATA 

- MOUXA 

♦OATA 

♦ OATA 

♦ OATA 

♦ OATA 

♦ CATA 
♦DATA 

♦DATA 

♦ OATA 

♦ OATA 

♦ OATA 
♦OATA 

•JgATjA 

♦ DATA 

♦ CATA 

♦ OATA 
♦OATA 

♦ DATA 
40 A?. A 

♦ DATA 

♦ CATA 

♦ DATA 

♦ OATA 

♦ DATA 
♦fiATA. 

♦ DA TA 

♦ DATA 

♦ OATA 

♦ OATA 
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IFPOT CIRCS RE/C 

CATA #h tl-fc SF AC (CRAFT 6 CCMRQi. S*£Ii_M _ 

c / f a • icne o itafe o 

C/T/4 IhrtEF l 

CAT A# IFCACC i 

C IT A* 1KPLR l 

C/TA# AS 7036.962 £ Q.04«9 E j 69.989 A *75.456 B» 160.798 F 270.0 
CATA* AS 711c. £ O.OCS El 9C.0 W'TJ.O 8% OVO F 0.0 

C/TA* I C AT 6 6IOei7 

C AT A ♦ IFtCT 1 _ . . _ . 

CATA# TIKE 0.0 1ST UP 20.0 FFEO 2.0 
CATA# CELTAT 0 .0 1 

C/TA# KFLCTS<1*9> 1 — - 

LATA# XFLCTS411) 101 OOOIOIOOO 1 0 1 0 0 0 I C l 000101000 

U AT A # KPLCT S ( 3 5 ) 1010CC101 000 

CATA# KPLCTS (41) 1 0 1 0 C C 

CATA# 1CR/V l 

CATA# PCCCS 1 1 l l 0 0 

CATA# MCDES l l l 1 1 1 0 

CATA# MOES 0 O O 0 0 0 

: at a# term tt .75 o.o o.o o.c o.o < u o — 

CATA* SCC 26*0 
CATA# 2*5 14.0 

l AT A* NELMS 7 

CATA# 2 LO it .0 it.O 36. C 36.0 C. 1 0.1 ?0.0 1.47 

CATA# 2LA 0.0 0.0 0,0 O.C 0 .C Q.O 

CATA# Z(2(l.l)2 1 . A 1 4 -1.414 -1.414 1.414 -1.6 -1.6 -2.0 

CATA# 262(2.1)2 0.1 0.1 0.1 0.1 2.0 -2.0 0.0 

CATA# 26 2(3. 1)2 - 1.4)4 -1.414 1.414 1.414 0.0 0.0 0.54 

C ATA# ALT AEK ( 1-7 ) O .0 

i AT A* FET/EK 45.0 13* .0 2 2 5.0 -45. C 0.0 0.0 180.0 
CATA* CANAEK 0.0 0.0 0.0 C.O 90.0 -90.0 -5.5 

CATA# E»COLSU-c) l.lStC*C7 
CATA# E*CCLS(l-t) l.4130*C7 
CATA# EROCLSl?) I.34P0*C7 
C ATA# FTUEF ' 1-C ) 0.5(2 

CAT A* FTU E E ( 7) 4.5 

CATA# FTUEEI1-6) O.OCOj 
CAT A# FTUEE (7) 0.QC2 

CATA# TFEFMtl-6) 8.t£D-0( 

;at\# TiP#s(j-e) i • <« 2 ou - o a 
C ATA* T IF #S ( 7 ) 0.2 6< 

JATA# SAC ( 1—7 ) 0.094 
CATA# STNK ( 1-6 > C.O 

CATA* STAK (7 ) 0*44 

CATA# Ff-CK ( 1-C ) 0 .002 l 
C/TA# PFCM 7) 0.C074 

CATA* FCC 0.0 
CATA* FCC 1.89C-07 
3 AT A* AEFt 0.0 
CATA* AC f C 2*2 

CATA? 1CAAF O IVlSCS 1 1ATTDF 0 IAFCCL 0 IMGMTS 0 

CATA* ICAfcP 0 1V1SCS l 1A X 4-CC- C 14# € 0- 1 IM G MS O 

CATA* 1 C A* P 0 IV ISCb 0 IATTC6 C InHEtL 1 IMGHTS 0 

OAT A# SK Q A (5.1 ) l 0.0 O.C 

DATA# CCA»PU.1>2 0.005 O.CC* C. CCS C. 0C5 0.005 0.005 

CATA# UFNC* 1 
CAT#* 16El*C# 0 

CATA# Ml.lll 0.0 0*0 Q.C 0.3 O.C 0.0 

C A T « # 2ETIAC 0.0 

CATA* ETT At O .0 - 

JATA* ETTAC 90.0 
CATA* UENSTY 5t.2 
CAT*# FACIEI 0.4(5 
ZATA* FAOFKG 1C.E2 
CATA* VISC^Y 2.0 
CATA* X#C ► IN (3 ) 0 .0 

CATA* XML * IM 1 ) 0.0 

CATA* CTCt 1.0 

CAT#* CMC 0.0 O.C 0.0 

CATA* I A C FT 2 ALFAt 0.0 EETAE -60.0 GA4AL 0.0 

OAIA* INCFT 2 ALFA6 C.O 861*8 0.0 <3 A * A f5 0.0" 

CATA* INQFT 2 ALF AE 0.0 3E T A r : 1 €0 . 0 GAMAF 190.0 

CATA** 1NTECFATUF MfcSSAGES 

CATA* 1MEEL 1 XMCMlN(2) .53C5 VMU 4500. DVM0M(21 0.0 

CATA** KNTPL(l) 0--NCTFING 1 — 2H0 0 ROEH SENSORS 2 — 4TH CPOFP 

CATA** KN11U2J 0 — PITCH CALX 

CATA** KNTFH3) 0 — NONE l — TACHOMETER 1ST CRDER DYNAMICS 
CAT** KSTFL 2 0 1 

CATA* KNT FL (9 ) 135791 — 

CATA# KMFL 1 0 1 

CATA* KM FL 1 0 C 

C/TA# KNTFL(IO) 3 ~ 


♦ CAT A 

♦ CATA 

♦ OAT A 

♦ C AT A 

♦ CATA 

♦ C AT A 

♦ * AT A 

♦ OAT A 

♦ OAT A 

♦ CAT A 

♦ CATA 

♦ CATA 

♦ CAT A 

♦ CAT A 

♦ CAT A 

♦ CAT * 

♦ CATA 

♦ CATA 

♦ CAT A 

♦ CAT A 

♦ CAT A 

♦ OAT A 

♦ CATA 

♦ CAT 4 

♦ OAT A 

♦ DAT A 

♦ CAT A 

♦ OAT A 

♦ CAT A 

♦ CAT A 

♦ OAT 4 

♦ OAT A 

♦ CAT A 

♦ OAT A 

♦ CATA 

♦ CAT A 

♦ CAT A 

♦ OAT A 

♦ D AT A 

♦ OAT A 

♦ CATA 

♦ CAT A 

♦ CAT A 

♦ OAT A 

♦ CAT A 

♦ DAT A 

♦ CAT A 
< CAT A 

♦ CAT A 

♦ DAT A 

♦ DAT A 

♦ OAT A 

♦ CAT A 

♦ CAT A 

♦ C*T 4 

♦ 0 4T A 

♦ OAT fi 

♦ OAT A 

♦ D M A 

♦ CAT 4 

♦ CAT * 

♦ CAT A 

♦ DAT A 

♦ OAT A 

♦ OAT A 

♦ CAT A 

♦ CAT A 

♦ CAT A 

♦ DAT A 

♦ DAT A 

♦ CAT A 

♦ CAT A 

♦ CAT A 

♦ DAT A 

♦ DAT A 

♦ r at a 

♦ DA T A 

♦ CAF A 

♦ CAT A 

♦ C/TA 

♦ OAT A 

♦ DAT A 
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OF POOR QUALi t V 


DE-B SPACECRAFT & CONTROL SYSTEM INPUT (CONTINUED) 


CJ1A+ KKTPl(lO) C 40ATA 

CATA4 CP«M(J|) 0.0* O.Q2- ♦DATA 

CAT A* (M((l) 0*0 ♦DATA 

2ATA4H FCCPTH CftCtH SENSlTH » I Th T ACHGME TEH 4 D AT A 

CAT A* (FAR* I*b tC.w. 2£4JO^. * I- ■!« — U 4 U SA. E - .000 3 ? . 1 6 * 0 6 OS 0* , 669 l 7 .1-q»X A CAT A 

CATA4 (FAPM15I t.CS Q * 0 A 1 6 0*0 0*4 0*0 C.O 1*0 * C AT A 

CliAA SVCS(I) -1.04? 0.0 -1.C4? C.O * 0 AT A 

C AT A* SVCS <1-20 1 0.0 Tilt -TTTT7 1 1 * > — HT^5 * CAT A 

CATA4 SV<S<|-20J 0.0 (111 C.O (If) 7C.5 * C AT A 

CATA* CP ARM22) Q .0349 ♦ CAT A 

CATA6 (PAt»t £2) C.2C3 -O.CC? 4CATA 

CAT*4 KEL(TS<2 16 J I 0 C 0 . 1 4 D AT A 

CAT A 4 KFLCTS(2£6J 1 C C 0-1 — — *CAT A 

CATA4 KFtCTS ( £ |6-i 1*0 ) 0 4 OAT A 

CATA4 CPAPN ( £ 1 1 -1*0 ACATA 

CATA4 <PAPtt(<l> 1*0 — ' 4CATA 

CATA4 <WCS (11) 1.5 A O AT A 

C AT A 4 S*CS(ll) O.C * DAT A 

CATA4 SVCSC12) 6 • c9 4 DAT A 

C IT A 4 EUFCCA(l) 0*0 0.0 O.C 4 CAT A 

- 4 1 A4 HUECtAH) JO.O £c.4 2 0* . O 4 CAT A 

C A T A 4 EUt<P(l) 0.0 -0.9 O.C 4 OAT A 

C4TA4 HUECP(l) 0.0 -1.2 0.;* 4 oat A 

C AT A 4 I CP AC 0 — — *OATA 

CATA4 ICP*C l 4 OAT A 

C AT A 4 ICPAC 2 4 OAT A 

CATA4 T C ( 1 ) 103S.778 IU27 # 0<5 iC^I.Cf* 5 4 C A'» A 

CAT A* KP < 1*1)2 20 20 20 20 £C 20 20 20 4C~TA 

CA1A4 CPA PA (41) -2 • l** 4 CAT A 

CATA4 KPLCT Sit 10 ) 000 4 DAT A 

C AT A 4 I AC AT S 0 4 DAT A 

CATA4 CPKACC2) 14500C.0 4 OAT A 

C A f A4 CCftS(l) 145000.0 4 DAT A 

C AT A 4 KELCTS(20QJ l 4 CAT A 

CA1A4 KPLCTS (2 16) 1 4 0 AT A 

l A T A 4 SVCS(l) 2.14 4 DAT A 

CAIA4 LCPU 40 4 DAT A 

C AT A 4 HL AST 0 4 OAT A 

C AT A * 1 4 DAT A 


vehsick apr isei 
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